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PREFACE 


This  volume  contains  all  the  Information  that  Is  required  to 
operate  the  computer  simulation  model  SNOHO.  The  Input  files  required, 
the  operational  logic  of  SNOMO  and  the  output  produced  by  SNOMO  are 
discussed  in  detail.  A  full  listing  of  the  SNOMO  Fortran-77  code  is 
Included. 
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CHAPTER  1 :  INPUT 


1.1  Cell  division. 


Acce'Sion  Fof 


NHS  CRAAI 
DIIC  TAb 

U-idf'io;'  'i.L-d 

JilSflf 


□ 

□ 


By 

Oibt'ibi.iioi'  / 


AvijildDility  Codes 


Dist 


Avdil  tnd  /  Of 
SpfCldl 


The  catchment  is  subdivided  into  computational  units,  called 
cells.  The  factors  used  for  this  subdivision  are  slope  angle, 
aspect,  elevation  and  vegetation  cover  type.  The  method  for 
catchment  subdivision  is  described  in  detail  in  chapter  5  (volume  1) 
and  need  not  be  repeated  here.  The  catchment  to  be  modelled  will 
therefore  consist  of,  after  subdivision,  a  number  of  cells  each  with 
an  aspect,  slope  angle,  elevation  and  vegetation  cover  value.  The 
location  of  each  cell  in  the  catchment  is  also  known.  The  depth  and 
occurrence  of  snowcover  is  calculated  for  each  cell  and  the  results 
for  all  the  cells  are  examined  together  in  order  to  calculate  the 
snowcover  depth  and  distribution  over  the  whole  catchment. 


1.2  Input  files. 

The  snowcover  calculations  are  conducted  separately  for  each 
cell.  The  presence  of  five  data  files  (snomo.dat ,  conLf .dat , 
decid.dat,  mix. dat  and  lapse.dat)  is  required  in  order  for  SNOHO  to 
operate  correctly.  The  input  required  for  these  5  files  will  be 
considered  in  detail  below  and  has  also  been  discussed  in  chapter  4 
(volume  1). 

1.2.1  snomo.dat. 

Tables  1.1,  1.2  and  1.3  show  sample  snomo.dat  files  for  cells  1 
and  24  in  1988  and  cell  1  in  1989  for  the  V3  catchment,  Sleepers 
River  Research  Watershed ,  Danville,  Vermont.  Table  1.4  shows  the 
variables  contained  in  snomo.dat .  The  variables  are  considered  in 
detail  below  in  the  order  that  they  appear  In  the  data  file. 
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Table  1.1.  Sample  snoioo.dat  file  for  cell  1,  U3,  1988. 


1  1  1988 

2  0.0  0.06  0.03  0.09  0.93  0.136 

3  0.0  0.24  0.08  0.94  0.40  0.410 

4  0.0  0.36  0.09  0.92  0.39  1.600 

3  1010  3  100 

6  6.13  -92.83  64  44.47 

7  1 

8  3 


9 

69.0 

10 

3.0 

11 

100.0  0.0  -1.0 

12 

39 

13 

64 

1100  -11.1 

-2.8 

71.00 

0.5 

0.7 

0.0 

14 

63 

1100  -16.3 

-3.8 

36.00 

0.3 

4.2 

0.2 

13 

66 

1100  -16.6 

0.0 

31.00 

0.3 

1.7 

0.0 

16 

67 

1100  -6.0 

3.6 

46.00 

0.7 

2.3 

1.3 

17 

68 

1100  -14.0 

4.4 

32.00 

0.0 

1.7 

0.1 

18 

69 

1100  -1.9 

4.6 

33.00 

1.0 

1.8 

6.1 

19 

70 

1100  -7.6 

2.4 

75.00 

0.3 

3.1 

1.6 

20 

71 

1100  -13.3 

-7.8 

42.00 

0.0 

3.4 

1.0 

21 

72 

1100  -16.9 

-1.6 

37.00 

0.5 

2.0 

0.0 

22 

73 

1100  -3.6 

3.6 

80.00 

0.3 

1.1 

4.6 

23 

74 

1100  -3.9 

0.3 

77.00 

1.0 

1.9 

1.4 

24 

73 

1100  -6.2 

-2.1 

73.00 

0.3 

2.2 

1.2 

23 

76 

1100  -6.3 

0.0 

68.00 

0.7 

1.9 

2.0 

26 

77 

1100  -10.9 

2.7 

36.00 

0.3 

2.3 

0.1 

27 

76 

1100  -9.7 

2.8 

43.00 

1.0 

1.2 

0.0 

28 

79 

1100  -13.3 

0.3 

46.00 

0.3 

2.8 

0.5 

29 

80 

1100  -19.0 

-6.1 

67.00 

0.3 

2.3 

2.8 

30 

81 

1100  -22.7 

-14.9  31.00  0.4  3.6  0.5 

31 

82 

1100  -19.1 

-*•2 

46.00 

0.3 

2.6 

0.0 

32 

83 

1100  -13.4 

4.8 

32.00 

0.4 

2.7 

0.0 

33 

84 

1100  -0.3 

11.0 

63.00 

0.3 

1.2 

0.0 

34 

83 

1100  -0.4 

7.2 

61.00 

0.3 

1.4 

0.0 

33 

86 

1100  4.4 

7.9 

92.00 

0.3 

0.8 

20.7 

36 

87 

1100  -2.2 

4.6 

92.00 

0.3 

3.5 

10.7 

37 

88 

1100  -6.1 

-0.2 

78.00 

0.3 

2.1 

1.2 

38 

89 

1100  -6.0 

12.7 

43.00 

0.3 

1.1 

0.0 

39 

90 

1100  1.8 

13.5 

31.00 

0.0 

2.1 

0.0 

40 

91 

1100  0.3 

10.8 

60.00 

0.0 

2.1 

0.0 

41 

92 

1100  0.3 

11.0 

38.00 

0.3 

2.2 

0.0 

42 

93 

1100  1.1 

16.8 

47.00 

0.3 

1.4 

0.0 

43 

94 

1100  4.3 

6.3 

88.00 

0.3 

1.4 

0.2 

44 

93 

1100  3.3 

11.8 

93.00 

1.0 

0.3 

0.4 

43 

96 

1100  3.3 

9.4 

93.00 

1.0 

1.4 

16.3 

46 

97 

1100  0.9 

13.3 

74.00 

0.6 

1.0 

0.0 

continuad  ovarlaaf . . . 
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47 

98 

1100 

0.3 

12.9 

61.00 

0.3 

1.0 

0.0 

48 

99 

1100 

3.4 

8.9 

87.00 

1.0 

1.1 

0.3 

49 

100 

1100 

1.0 

8.7 

43.00 

0.3 

2.9 

0.3 

SO 

101 

1100 

-0.3 

8.7 

82.00 

0.0 

3.2 

0.0 

31 

102 

1100 

-0.7 

8.0 

76.00 

0.3 

1.8 

1.2 

32 

103 

1100 

-0.3 

12.3 

33.00 

0.3 

1.2 

0.0 

Interpretation  of  snomo.dat  file; 

LINE  1  call  Idantlfieation  niaibaE,  simulation  yaar. 

2  -1  -2  -1  -1 

LINE  2  0.0,  thanaal  dlffuslvlty  (en  ailn  ).  hast  conductivity  (calcm  min  ‘K  ),  aoissivity 

-3 

(daciaial),  albado  (daciiaal),  danalty  (gm  );  lor  naw  snow. 

LINE  3  as  lina  2,  but  lor  old  snow. 

LINE  4  as  lina  2  and  3,  but  for  sandy  soil. 

LINE  S  air  prassura  (ob),  cloud  typa,  instruosnt  haight  (cb). 

LINE  6  slopa  angla,  aspaet,  start  data,  Istituda. 

LINE  7  lapsa  rata  indicator. 

LINE  8  vagatatlon  eovar  typa  indicator. 

LINE  9  Initial  anowdapth,  (cai). 

LINE  10  critical  anowdapth  (5em). 

LINE  11  dapth  of  aoil  protila  (ca),  soil  tsoparatura  at  dapth  (*C),  aurfaca  soil  tomperaturs 

<*C). 

LINE  12  niaabar  of  days  in  data  fila  Binua  ana. 

LINES  13  daily  SMtaorological  data:  Julian  data,  obsarvation  tiow,  atlnlBUB  air  taaparatura 
(*C),  Biaaiaaaa  air  taoparatura  (*C), 

TO  32  ralativa  bunidity  (X),  cloud  cover  (0-1),  wind  spsad  (bs  ^). 


Tabla  1.2 


Sample  snoao.dat  file  for  cell  24 


W3 


1988. 


1  24  1988 

2  0.0  0.06  0.03  0.09  0.93  0.136 

3  0.0  0.24  0.08  0.94  0.40  0.410 

4  0.0  0.36  0.09  0.92  0.39  1.600 
3  1010  3  100 

6  6.27  -101.37  64  44.47 

7  0 

8  0 

9  69.0 

10  3.0 

11  100.0  0.0  -1.0 

12  39 

13  64  1100  -11.1  -2.8  71.00  0.3  0.7  0.0 
Linas  14  to  32  as  flla  (1). 


Table  1.3.  Sample  snomo.dat  file  for  cell  1,  U3,  1989 


1  1  I860 

2  0.0  0.06  0.03  0.98  0.9S  0.136 

3  0.0  0.24  0.06  0.84  0.40  0.410 

4  0.0  0.36  0.09  0.92  0.39  1.600 
3  1010  3  100 

6  6.13  -92.63  64  44.47 

7  1 

8  3 

9  32.9 

10  3.0 

11  100.0  0.0  -1.0 

12  46 


13 

64 

1100  -8.3 

0.7 

96.00 

0.3 

0.3 

7.1 

14 

63 

1100  -19.9 

-9.4 

61.00 

1.0 

2.3 

0.0 

13 

66 

1100  -28.7 

-13.8 

63.00 

0.3 

1.7 

0.0 

16 

67 

1100  -20.4 

-3.9 

72.00 

0.0 

1.3 

0.0 

17 

68 

1100  -14.9 

-1.2 

67.00 

0.0 

1.4 

0.0 

18 

69 

1100  -11.3 

1.6 

32.00 

0.3 

1.7 

0.0 

19 

70 

1100  -7.6 

2.02 

47.00 

0.3 

1.9 

1.02 

20 

71 

1100  -12.8 

-1.7 

87.00 

0.3 

2.4 

1.6 

21 

72 

1100  -16.0 

-1.3 

40.00 

0.0 

1.2 

0.0 

22 

73 

1100  -3.9 

0.6 

61.00 

0.3 

1.0 

0.0 

23 

74 

1100  3.4 

14.3 

96.00 

0.3 

1.8 

0.0 

24 

73 

1100  -4.6 

2.3 

78.00 

0.3 

2.2 

0.0 

23 

76 

1100  -6.7 

2.6 

73.00 

1.0 

2.1 

3.1 

26 

77 

1100  -8.1 

1.0 

97.00 

0.3 

0.4 

17.0 

27 

78 

1100  -12.3 

-6.1 

64.09 

0.3 

1.7 

1.02 

28 

79 

1100  -12.1 

-0.8 

87.00 

1.0 

1.2 

1.32 

29 

80 

1100  -9.3 

-0.9 

94.00 

0.3 

2.0 

7.6 

30 

81 

1100  -14.3 

-3.6 

37.00 

0.3 

1.8 

0.0 

31 

82 

1100  -10.7 

2.6 

47.00 

0.0 

1.8 

0.0 

32 

63 

1100  -9.4 

2.2 

53.00 

0.3 

1.8 

2.3 

33 

84 

1100  -2.7 

6.3 

68.00 

0.3 

0.8 

3.1 

34 

83 

1100  0.0 

6.3 

91.00 

0.3 

1.1 

0.0 

33 

86 

1100  0.0 

16.9 

60.00 

0.3 

0.9 

0.0 

36 

87 

1100  7.0 

17.1 

79.00 

0.2 

0.8 

2.3 

37 

88 

1100  -0.2 

13.6 

63.00 

1.0 

1.7 

13.3 

38 

89 

1100  -2.4 

0.4 

93.00 

0.3 

1.2 

18.3 

39 

90 

1100  -2.7 

0.3 

97.00 

1.0 

0.6 

10.2 

40 

91 

1100  -2.3 

0.2 

94.00 

0.3 

1.8 

4.3 

41 

92 

1100  -3.4 

6.3 

83.00 

0.3 

1.3 

0.0 

42 

93 

1100  1.0 

3.6 

94.00 

1.0 

1.1 

26.2 

43 

94 

1100  1.2 

4.6 

99.00 

1.0 

0.8 

6.9 

44 

95 

1100  4.3 

14.0 

90.00 

0.4 

0.9 

8.4 

43 

96 

1100  0.1 

6.3 

96.00 

0.3 

1.8  21.3 

46 

97 

1100  0.1 

4.2 

94.00 

0.6 

1.7 

9.4 

eontlnuad  ov»rl«af . . 
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47 

98  1100  '3.4 

2.6 

82.00 

0.3 

2.4 

0.0 

48 

99  1100  -4.9 

6.2 

43.00 

0.3 

1.3 

3.1 

49 

100 

1100  -S.7 

1.2 

70.00 

0.3 

1.6 

6.6 

SO 

101 

1100  -7.9 

-0.2 

70.00 

0.7 

1.7 

0.0 

31 

102 

O 

O 

1 

4.2 

66.00 

0.3 

0.7 

0.0 

S2 

103 

1100  -1.1 

S.9 

S8.00 

0.3 

1.9 

0.0 

S3 

104 

1100  -0.3 

S.2 

73.00 

0.9 

1.1 

0.0 

S4 

lOS 

1100  0.4 

6.8 

71.00 

0.3 

2.1 

6.4 

SS 

106 

1100  1.0 

4.6 

97.00 

0.3 

0.9 

4.6 

S6 

107 

1100  3.S 

14.7 

80.00 

0.3 

0.8 

0.0 

S7 

108 

1100  0.4 

11.1 

93.00 

0.3 

2.4 

3.8 

S8 

109 

1100  -0.9 

S.9 

33.00 

0.3 

1.9 

0.0 

S9 

110 

1100  -0.6 

3.2 

90.00 

0.3 

2.1 

1.8 

7 


Table 


LIME 


1 


2 


3 


4 


S 


6 


7 


S 


1.4.  Description  of  variables  in  snomo.dat. 


VARIABU  SXAL/ 
IHTBOBR 


DBSCRIPTIOe 


CELL 

I 

M 

• 

! 

of  call. 

YEAR 

I 

Year  of  aiaulation 

VHEHSIld.l) 

R 

Hom  snow  vector,  0.0. 

VIIEUSH(1.2) 

R 

- 

**  ,  tbttmal  dlffusivlty  (cm  min  ^ 

VKEWSMd.S) 

R 

*• 

"  .  heat  conductivity  (calca  ^ain 

VHEMSMd.S) 

R 

s. 

"  .  aaissivity  (daclaal). 

VHEHSHd.S) 

R 

- 

"  ,  albedo  (declael). 

VHEHSRd.6) 

R 

**  »* 

"  ,  density  (gca  ^). 

VQLOSNd.l) 

R 

Old  snow  vector,  0.0. 

VOLDSHd.2) 

R 

*s  as 

**  ,  thnimal  dlffusivlty  (cm  mln’^l 

VOLDSHd.3) 

R 

•*  " 

"  ,  heat  conductivity  (calca  *Bln' 

VOLDSHd.i) 

R 

~  ,  aaissivity  (deciaal). 

TOLDSHd.S) 

R 

"  ,  albedo  (deciaal) . 

VOLDSMd.g) 

R 

«  ** 

,  density  (gca  ^). 

VSOILd.l) 

R 

Soil  vector,  0.0. 

VSOILd.2) 

R 

"  ,  thszmsi  dlffusivlty  (ca^sdn 

VSOILd,3) 

R 

"  ,  heat  conductivity  (calca  ^ain  ^’K 

VSOILd.l) 

R 

,  aaissivity  (deciaal) . 

VSOILCl.S) 

R 

*• 

"  ,  albedo  (daclaal). 

VSOILd.S) 

R 

** 

-3 

"  ,  density  (gca  ) . 

PRESS 

I 

Air  pressure  (d>). 

HCLOUO 

I 

Cloud  typo,  see  table  1.2. 

ZA 

I 

Instrunent  height  (ca). 

SLOPEl 

R 

Call  slope  angle  (*) 

stntFCi 

R 

Coll  aspect  (*) 

DAYl 

I 

Start 

date  (Julian  calender). 

LAY 

R 

Latitude  of  catefanant  (*). 

ILAFSE 


rVEG 

rVEa~l.  conifarous  etmz . 

IVC0~2,  deciduous  cover. 

IVE0“3,  Klzed  cover. 

tVEO^O,  pasture  sad  cLaar-cut  (open). 


I  Lapse  rata  uulicator. 

ILAPSE~0,  elevatloo  band  ISOO- 2000ft. 
ILAPSE'l,  elevation  band  lOOO-lSOOft. 

t  Call  vegetation  cover  tppa  indicator. 


continued  overleaf. . . 
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9  SIIDFl  K 

10  CRI^'  It 


Initial  anovdapth  (ca) . 
Critical  anowdapth  (Sea). 


11  SOILOP  R 

SOILTP  R 

SSOLTP  R 


Oapth  of  aoil  profila  (ea). 

Soil  taaparatura  at  dapth  (*C>. 
Suzfaea  aoil  taiqMratura  (*C). 


12  M 


I  HurtMz  of  daya  to  bo  aadallad.  Muabat  of  daily 

taatoorolosical  data  daya  alnua  1. 


XXXd.a)  I 

md.l)  I 

YYYd.l)  R 

YYYd.9)  R 

mfd,2)  R 

YYYd,3)  1 

md,6)  R 

md.7>  R 


Dally  aMtaorolotlcal  data. 

Data  (Julian  ealandar). 

Qbaarvatlon  tiaa  (24hr  clock). 
Hiniaaa  daily  air  taaparatura  (‘O. 
Maziaiaa  dally  air  taaparatura  (‘O. 
Ralativa  huaaldlty  (X). 

Cloud  cover  (0-1). 

Wind  apaad  (aa  ^). 

Praeipitation  (b  watar). 
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LI££_L 

1.  CELL 

The  number  of  the  cell.  Each  cell  is  numbered  for  easy 
identification. 

2.  YEAR 

The  year  of  the  simulation. 

LINE  2:  'New'  snow  vector. 

The  vector  for  the  physical  characteristics  of  'new'  snow. 

Chapter  4  (volume  1)  considered  the  division  of  the  physical 

characteristics  of  snow  into  those  commonly  associated  with  'new' 

(fresh)  snow  and  'old'  snow.  Values  are  taken  from  various  sources 

or  measurements  on  site  if  available.  The  physical  characteristics 

required  are  thermal  dlffusivity  (cm  min  ) ,  heat  conductivity 
•  2  - 1  1 

(calcm  min  K  ),  emissivity  (decimal),  albedo  (decimal)  and 

.  3 

density  (gem  ) .  The  first  value  in  the  vector  is  always  input  as 
zero.  This  is  because  it  represents  snowpack  depth  and  is  calculated 
within  the  program. 

LINE  3:  'Old'  snow  vector. 

As  line  2,  except  concerned  with  the  values  of  the  physical 
characteristics  of  'old'  snow. 

LINE  4:  Soil  vector. 

As  lines  2  and  3,  except  concerned  with  the  values  of  the 
physical  characteristics  of  the  soil.  Soil  type  is  taken  as  the 
predominant  soil  type  for  that  particular  cell  or  catchment.  The 
soil  vector  is  rarely  used  and  therefore  the  accuracy  of  the  soil 
values  is  not  particularly  important.  Again  values  can  be  taken  from 
the  literature,  or  on  site. 
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LI££L1 

1 .  PRESS . 

Air  pressure  (mb).  The  average  air  pressure  for  the  simulation 
period  is  used,  or  a  median  value. 

2.  NCLOUD. 

Values  for  cloud  cover  type  are  contained  in  table  1.5.  The 
average  cloud  cover  type  is  used  or,  in  the  absence  of  sufficient 
data,  a  default  cloud  cover  type  of  5  (stratocumulus)  is  used. 

3.  ZA. 

Instrument  height  is  the  average  height  (cm)  of  the  instruments 
used  to  measure  the  meteorological  data,  ie.  wind  speed  and  relative 
humidity,  from  the  grotmd/snow  surface.  If  this  measurement  is 
unavailable  or  is  too  variable  a  default  value  of  100cm  is  used. 

LIWE-§ 

1.  SLOPEl. 

Cell  slope  angle  (*). 

2.  SURFCl. 

Cell  aspect.  The  cell  aspect  is  measured  in  degrees  from  South 
with  westerly  values  positive  and  easterly  values  negative,  ie. 
Hest-+90*  and  East— 90*. 

3 .  DAYl . 

Date  of  start  of  simulation,  Julian  calender. 

4.  LAT. 

Latitude  of  catchment  or  cell  (degrees  and  minutes) .  The 
latitude  of  the  midpoint  or  mouth  of  the  catchment  can  be  used,  at 
the  operator's  discretion. 

LUifLJL 

1.  ILAPSE. 

Lapse  rate  indicator.  SNOMO  is  currently  set  up  for  the  U3 
catchment.  Sleepers  River  Research  Watershed,  Danville,  Vermont. 

When  subdivided  the  elevations  derived  for  each  cell  ranged  within 
1000-2000ft  a.s.l..  A  midpoint  of  1500ft  a.s.l.  was  used  as  a 
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Table  1.5.  Cloud  Genera  and  Cloud  Type  Indices  ,  Balick  et.al. 
(1981a). 


Cloud  Abbrovlotlon  Index  CaoMnts 

Genera  value 

Cirrus  Ci  1  High  clouds  eoaposad  of  white  delicate  filasente,  patches 

of  narrow  bends,  alaDants  often  curved  or  slanted  and 
aaallar  than  Cs,  never  overcast  or  precipitatins. 


Cirrostratus  Ce  2  High  clouds  appearing  as  whitish  veil  usually  fibrous, 

often  produces  halo  phanooena,  thinner  than  As,  doss  not 
appear  to  Bx>ve,  nonprecipitating. 


Altocussilus  Ac  3  Midleval  clouds,  patches,  usually  broken,  lee  wave  clouds, 

elements  smeller  than  Sc,  nonprecipitating. 

Altostratus  As  4  Hldleval  gray  sheet  or  layer  of  striated,  fibrous  or 

uniform  appearance,  large  horisontal  extant;  thicker  than 
Cs,  thinner  than  Ha,  precipitation  generally  light  and 
continuous  (if  any). 


Stratoctxoulus  Sc 


Stratus  St 


Hlobostratus  Hs 


Fog  FG 


3  Grey  and/or  whitish  layer  or  patch,  nearly  always  has  dark 
spots  and  is  nonflbrous;  elsmants  larger  than  Ac, 
nonpraclpitatlng . 

6  Grey  rather  uniform  base,  patches  ragged  if  present, 
precipitation  imusual  but  light  and  continuous  if  present, 
lower  md  mors  uniform  than  Sc,  less  danse  and  less  ’wet’ 
than  Hs. 

7  Grey  often  dark,  diffuse,  large  horlxontal  and  vertical 
extant,  thicker  than  As,  more  uniform  than  Sc,  often 
precipitating,  precipitation  continuous. 

8 


*  Cloud  genera;  Cumulus  (Cu),  Clrrociaulus  (Cc)  and  Ciasulonlnbus  (Cb)  are  not  treated  hare. 
At  low  cloud  covers  (0.3)  Cu  and  Cc  may  be  approximated  with  Ac. 


k 
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further  subdivision  of  the  elevation  resulting  in  2  elevation  bands: 
1000-1500ft  and  1500-2000ft.  If  the  elevation  of  the  cell  is  between 
1000-1500ft  the  lapse  rate  indicator  is  set  to  1,  if  between  1500- 
2000  then  it  is  set  to  0.  The  elevation  bands  can  be  altered  if 
required. 

LINE  g 

1.  IVEG. 

Vegetation  cover  type  indicator.  Four  vegetation  cover  types 
can  be  modelled  using  SNOMO: 

IVEG— 1,  coniferous  cover 

IVEG-2,  deciduous  cover 

IVEG-3,  mixed  cover 

IVEG-0,  pasture  or  clear-cut  (open). 

The  vegetation  cover  types  are  those  found  at  U3.  Again,  as 
with  ILAPSE,  these  can  be  altered  if  necessary. 

LIKE  ? 

1.  SNDPl. 

Initial  snowdepth.  This  is  the  known  (measured  or  estimated) 
snowdepth  for  the  cell  at  the  start  of  the  simulation.  If  estimated 
the  value  for  initial  snowdepth  is  based  on  the  relationship  between 
snowdepth,  vegetation  cover  type  and  elevation.  If  one  value  at  a 
base  station  (1500-2000ft  a.s.l)  is  available  then: 

(1)  Cell  with  elevation  <  1500ft  a.s.l.,  -10cm. 

(2)  Cell  with  coniferous  forest  cover,  -5cm. 

(3)  Cell  with  mixed  forest  cover,  -t-5cm. 

(A)  Cell  with  deciduous  forest  cover,  +10cm. 

These  values  can  be  used  or  compiled  according  to  another  catchment. 

LlK£JLfi 

1.  CRITDP. 

Critical  depth  (5cm).  This  is  the  depth  below  which,  because 
of  the  difficulty  of  accurate  measurement,  no  snow  is  said  to  be 
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present. 

UPE  11 

1.  SOILDP. 

Depth  of  soil  profile  (cm).  Default  value  is  100cm. 

2.  SOILTP. 

Soil  temperature  at  depth  (*C).  Default  value  is  0*C.  SNOMO 

does  not  currently  model  the  temperatures  within  the  soil  or  the  heat 

flux  across  the  ground/snow  interface  (Q  )  and  therefore  the  accuracy 

c> 

of  SOILTP  is  not  vital. 

3.  SSOLTP. 

Surface  soil  temperature  (‘C).  Default  value  of  -1*C  is  used. 
Again,  for  the  same  reasons  as  SOILTP,  accuracy  is  not  vital. 

UHE  12 

1.  N 

The  number  of  days  to  be  modelled  in  the  simulation.  This  is 
calculated  as  the  number  of  daily  meteorological  data  available  in 
snomo.dat  minus  one. 

LINE  13  and  onward  :  Daily  meteorological  data. 

1.  XXX(1,8) 

Julian  date. 

2.  XXX(1,1) 

Observation  time  (24hr  clock),  default  value  of  llOOhrs. 

3.  YYY(1,1) 

Minimum  air  temperature  (*C). 

4.  YYY(1,9) 

Maximum  air  temperature  (*C). 

5.  YYY(1,2) 

Relative  humidity  (%). 

6.  YYY(1,3) 

Cloud  cover  amount  (decimal). 

7.  YYY(1,6) 
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Wind  speed  (ms‘^) 

8.  YYY(1,7) 

Precipitation  (mm  water) . 

Table  1.6  is  provided  for  the  guidance  of  the  compilation  of  a 
snomo.dat  file. 

1.2.2  conif.dat,  decid.dat  and  mix.dat. 

These  three  data  files  contain  the  information  required  if  the 
vegetation  cover  options  of  coniferous,  deciduous  or  mixed  cover  re 
used.  The  files  contain  values  of  the  variables  o^,  x,  Cf >  <*f  and 
Z£.  Examples  of  the  three  data  files  used  for  the  U3  simulation  are 
shown  in  table  1.7.  The  variables  contained  in  the  files  are 
considered  in  chapter  4  (volume  1)  and  below: 

1.  Of,  foliage  cover  fraction  (0-1). 

This  describes  the  density  of  the  vegetation  cover.  o^-O 
represents  no  foliage  and  therefore  no  radiative  shielding  and  (7£-l 
represents  complete  radiative  blocking.  Table  1.4  shows  various 
limiting  values  of  oj  taken  from  Geiger  (1965)  and  Deardorff  (1978) 
for  various  vegetation  covers. 

2.  x<  state  of  vegetation  (1-1000). 

X  is  used  as  a  multiplier  of  the  stomatal  resistance  function. 

A  summer  value,  when  the  vegetation  is  healthy  and  active,  is  1,  a 
winter  and  therefore  dominant  or  dead  vegetation  cover  value  is  1000. 
Other  values  can  be  chosen  to  adjust  stomatal  resistance  for  moisture 
stress,  senescence  or  other  factors. 

3.  <£,  foliage  emisslvity  (decimal). 

4.  a£,  foliage  albedo  (decimal). 

5.  Z£,  foliage  hei^t  (cm). 

Values  for  <£  and  can  be  obtained  from  the  literature. 
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Table  1.6.  Compilation  guideline  for  snoiao.dat. 


All  input  on  the  same  line  is  separated  by  at  least  1  space  and 
starts  in  the  first  column. 


LINE  1 

1.  Number  of  cell 

2.  Year  of  simulation 

LINE  2 

1.  0.0 

2.  New  snow  thermal  diffusivity 

3.  "  "  heat  conductivity 

4.  "  "  emissivlty 

5 .  "  "  albedo 

6.  "  "  density 

LINE  3 

1.  0.0 

2.  Old  snow  thermal  diffusivity 

3.  "  "  heat  conductivity 

4.  "  "  emissivlty 

5.  "  "  albedo 

6.  "  "  density 

LINE  4 

1.  0.0 

2.  Soil  thermal  diffusivity 

3.  "  heat  conductivity 

4 .  "  emissivlty 

5 .  "  albedo 

6.  "  density 


no 


continued  overleaf . . . 


LINE  5 

1 .  Air  pressure 

2 .  Cloud  type 

3.  Instrument  height 

LINE  6 

1.  Cell  slope  angle 

2.  Cell  aspect 

3.  Start  date 

4.  Latitude  of  catchment  or  cell 

LINE  7 

1.  Lapse  rate  indicator 

LINE  8 

1.  Cell  vegetation  cover  type  indicator 
LINE  9 

1.  Initial  snowdepth 

LINE  10 

1.  Critical  snowdepth 

LINE  11 

1.  Depth  of  soil  profile 

2.  Soil  temperature  at  depth 

3.  Surface  soil  temperature 

LINE  12 

1.  Number  of  days  to  be  modelled 


continued  overleaf. . . 


1.  Julian  date 

2 .  Observation  time 

3.  Minimum  daily  air  temperature 

4.  Maximum  daily  air  temperature 

5.  Relative  humidity 

6.  Cloud  cover  amount 

7 .  Vind  speed 

6,  Precipitation 


Table  1.7.  conif.dat,  decid.dst  and  mix.dat  files  used  for  U3 

conLf  .dat\  q. 7a,  iooo.o.m, 0.10,2000 

decid.dat:  0.  so,  1000, 0.07,0.  is, 2000 

mix.dat:  0.60,1000, 0.97s, 0.12s, 2000 

Interpretation: 

1.  0£,  foliage  cover  fraction  (0-1) 

2.  Xi  state  of  vegetation  (1-1000) 

3.  <£,  foliage  emissivity  (decimal) 

4.  a£,  foliage  albedo  (decimal) 

5.  Z£,  foliage  height  (cm). 


Table  1.8,  lapse.dat  file  used  for  V3. 
lapse. dat:  2. 0,1.0 

Interpretation: 

1.  TMAX,  lapse  rate  factor  for  dally  maximum  temperature  CC) 

2.  THIN,  lapse  rate  factor  for  daily  minimum  temperature  CC) 
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1.2.3  lapoe-dat. 

The  file  lapse.dat  contains  values  for  TMAX  and  TMIN.  These 
are  the  lapse  rate  alteration  of  the  maximum  (TMAX)  and  minimum 
(TMIN)  dally  air  temperatures.  The  values  for  the  lapse  rate 
alteration  are  those  used  for  the  catchment  U3  and  are  obtained  from 
the  literature  (chapters  4  and  6,  volxime  1).  If  applied  to  different 
catchments  the  values  can  be  altered  If  necessary.  An  example  of  the 
lapse.dat  file  used  for  the  U3  simulation  Is  shown  In  table  1.8. 
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CHAPTER  2:  OUTPUT 


2.1  Output  file. 

Table  2.1  shows  the  output  file  SN01K).RES  for  cell  1,  V3,  1988. 
The  input  file  used  was  that  shown  in  table  1.1.  The  output 
variables  are  considered  below  in  the  order  in  which  they  appear  in 
SNOMO.RES  and  using  the  titles  used  in  SNOMO.RES . 

1.  CELL  NO. 

The  identification  number  of  the  cell. 

2.  YEAR 

The  year  of  the  simulation. 

3.  JULIAN  DATE 

Date,  Julian  calender. 

4.  SNOUDEFTH 

Calculated  depth  of  snowpack  (cm) . 

5 .  SNOWMELT 

Calculated  dally  depth  of  snow  melted,  in  both  centimetres  of 
snow  and  millimetres  of  water  equivalent. 

6 .  SNOWFALL 

Dally  snowfall  (mm  water  equivalent) .  The  value  for  snowfall 
on  the  first  day  of  simulation  is  equal  to  the  initial  snowpack  depth 
and  is  not  related  to  the  amount  of  daily  snowfall. 

7 .  DENSITY 

Daily  average  snowpack  density  (gem  ). 

8.  RAIN- ON- SNOW 

Occurrence  of  a  rain-on- snow  event.  1-rain- on- snow  event,  0-no 


rain -on- snow  event. 
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Table  2.1.  Saaple  SNOMO.RES  file,  for  cell  1,  U3,  1988 


CZLZ.  no.:  1 
YEAR:  1988 


miAH 

SHOHOEITH 

SHCMIELT 

SMOUFAU. 

dehsity 

RAIH-OM-SHOH 

DATE 

(ea) 

(ca  snow) 

(M  watar) 

(Cffl) 

<G/H»*3) 

( l“yaa .  0”no ) 

64.00 

69.00 

0.00 

0.00 

60.00 

0.39 

0 

83.00 

69.13 

0.00 

0.00 

0.15 

0.39 

0 

68.00 

69.13 

0.00 

0.00 

0.00 

0.39 

0 

67.00 

66.81 

0.00 

0.00 

1.10 

0.39 

0 

68.00 

66.86 

0.00 

0.00 

0.07 

0.41 

0 

68.00 

71.37 

0.00 

0.00 

4.49 

0.41 

0 

70.00 

72.33 

0.00 

0.00 

1.18 

0.39 

0 

71.00 

73.28 

0.00 

0.00 

0.74 

0.38 

0 

72.00 

73.28 

0.00 

0.00 

0.00 

0.38 

0 

73.00 

71.73 

0.00 

0.00 

3.53 

0.38 

0 

74.00 

72.78 

0.00 

0.00 

1.03 

0.40 

0 

73.00 

73.66 

0.00 

0.00 

0.88 

0.39 

0 

76.00 

75.13 

0.00 

0.00 

1.47 

0.39 

0 

77.00 

73.29 

0.00 

0.00 

0.07 

0.38 

0 

78.00 

73.20 

0.00 

0.00 

0.00 

0.36 

0 

79.00 

70.80 

0.00 

0.00 

0.37 

0.38 

0 

80.00 

67.03 

3.91 

24.24 

2.06 

0.41 

0 

81.00 

61.97 

3.43 

22.34 

0.37 

0.41 

0 

82.00 

60.20 

1.77 

7.26 

0.00 

0.41 

0 

83.00 

60.20 

0.00 

0.00 

0.00 

0.41 

0 

84.00 

38.33 

1.83 

7.60 

0.00 

0.41 

0 

83.00 

33.32 

3.03 

12.42 

0.00 

0.41 

0 

86.00 

30.89 

4.43 

18.14 

0.00 

0.41 

1 

87.00 

34.10 

4.66 

19.13 

7.87 

0.41 

0 

88.00 

34.98 

0.00 

0.00 

0.88 

0.39 

0 

89.00 

34.98 

0.00 

0.00 

0.00 

0.39 

0 

90.00 

32.23 

0.00 

0.00 

0.00 

0.38 

0 

91.00 

47.38 

4.67 

19.14 

0.00 

0.41 

0 

92.00 

43.43 

4.16 

17.04 

0.00 

0.41 

0 

83.00 

38.68 

4.73 

19.46 

0.00 

0.41 

0 

94.00 

34.08 

4.60 

18.83 

0.00 

0.41 

1 

93.00 

30.26 

3.83 

13.69 

0.00 

0.41 

1 

96.00 

23.42 

4.84 

19.83 

0.00 

0.41 

1 

97.00 

19.96 

3.46 

22.38 

0.00 

0.41 

0 

98.00 

14.72 

3.24 

21.48 

0.00 

0.41 

0 

99.00 

9.87 

4.85 

19.88 

0.00 

0.41 

1 

100.00 

4.00 

3.88 

24.10 

0.00 

0.41 

1 
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CHAPTER  3:  PROGRAMME  DETAILS 


3.1.  Programme  structure. 

Chapter  A,  volume  1,  considers  In  detail  the  equations  used  by 
SNOMO  to  calculate  snowdepth  within  one  cell.  This  chapter  presents 
the  programme  logic  that  enables  these  equations  to  be  utilised  and 
to  Interact  successfully.  Figure  3.1  shows  the  SNOMO  programme 
logic.  SNOMO  broadly  consists  of  2  sections; 

(1)  Control  programme. 

This  Is  responsible  for  the  determination  of  he  type  of 
precipitation,  the  calculation  of  melt,  the  calculation  of  the  affect 
of  a  raln-on-snow  event,  the  calculation  of  the  heat  flux  across  the 
ground/snow  Interface  iQ  ) ,  various  operational  markers ,  the 
calculation  of  the  depth  of  the  snowpack  and  the  physical 
characteristics  of  the  snowpack,  the  model  Input  and  the  model 
output . 

(2)  Algorithm  to  calculate  the  majority  of  the  components  of  the 
snowpack  energy  budget. 

This  algorithm  Is  responsible  for  the  calculation  of  the 
snowpack  energy  budget  components  K1 ,  Kt ,  LI ,  LT ,  and  and  the 
snowpack  surface  and  Internal  temperatures.  The  algorithm  Is  adapted 
from  the  US  Corps  of  Engineers  Terrain  surface  Temperature  Model 
(TSTM)  as  developed  by  Ballck  et.al.  (1981a  &  b) .  The  basic 
equations  used  by  the  original  TSTM  remain  the  same  as  In  the 
modified  TSTM.  TSTM  has  been  heavily  modified  from  Its  original 
state  as  presented  in  Ballck  et.al.  (1981a  &  b)  and  has  also  been 
Incorporated  into  the  logic  structure  of  SNOMO.  The  logic  structure 
of  TSTM  is  shown  In  figure  3.2.  TSTM  is  operational  for  situations 
where  a  vegetation  cover  is  present  or  absent.  Vegetation  cover  Is 
modelled  using  a  vegetation  cover  algorithm,  VEGIE,  within  TSTM. 

The  relationship  between  TSTM  and  VEGIE  is  shown  in  figure  3.3. 


Slnpllfled  flowchart 


Figure  3.2.  Slatpllfied  flowchart  for  TSTM,  Balick  et.al. 
(sheet  2  of  S) . 


1981a 
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Figure  3. 


T 


2.  Sioplifled  flowchart  for  TSTM,  Ballck  et.al. .  1981a 
(sheet  3  of  5) . 
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Figure  3.3.  Sequence  of  calculation  of  the  major  components  of  the 
TSTM/VEGIE  system,  Balick  et.al.  (1981b) . 
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In  block  1  the  model  is  Incremented  one  time  step  or  Initialized,  or 
terminated)  in  accordance  with  procedures  established  for  TSTM.  In 
block  2  the  decision  %fhether  to  utilize  VEGIE  or  not  is  made.  If 
VEGIE  is  not  used  a  non-vegetated  surface  energy  budget  is  evaluated 
as  an  upper  boundary  condition  (block  3)  for  the  solution  of  the 
equation  of  heat  transfer  through  the  terrain  materials  (block  6) . 

The  solution  for  the  surface  temperature  (block  7)  comes  from  the 
evaluation  of  the  surface  energy  budget  equation  and  therefore  to 
solve  the  heat  transfer  equation.  In  order  to  achieve  this  the  heat 
conduction  term  and  the  distribution  of  heat  in  terrain  materials 
must  be  calculated  (block  6).  If  VEGIE  is  used,  the  energy  budget  of 
block  3  is  replaced  by  blocks  4  and  5.  Blocks  4  and  S  comprise 
VEGIE.  Block  4  calculates  the  energy  budget  for  the  foliage  that 
includes  a  contribution  from  the  ground  surface.  Block  5  calculates 
the  energy  budget  for  the  ground  surface  that  includes  a  contribution 
from  the  foliage  layer. 

Solutions  of  temperature  for  the  foliage  (from  block  4)  and  the 
ground  (from  block  5)  are  performed  by  a  simple  root-finding 
algorithm  and  are  combined  according  to  the  proportion  of  foliage 
cover  to  yield  an  average,  or  effective,  temperature  of  the  vegetated 
surface.  The  ground  energy  budget  is  then  used  in  the  evaluation  of 
heat  flow  in  the  terrain  (block  6),  and  the  programme  returns  to 
block  1. 

3.2.  Additional  information. 

SNOMO  can  be  altered  very  simply  in  order  to  produce,  in  the 
output  file,  the  hourly  energy  budget  values  for  the  snowpack  of  K4 , 
Kt,  LI,  Lt,  0^  and  Q^.  The  alterations  are  found  under  the  heading 
'OPTION  2'  in  the  hardcopy  of  the  programme  (section  3.3.) 

The  absence  of  a  mass  conservation  routine  was  discussed  in 
voltaie  1.  A  mass  conservation  routine  has  been  written  and  tested 
(chapter  7,  volume  1)  and  is  included  in  Che  hardcopy  of  the 
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prograooe  (section  3.3.)  as  'OPTION  3'. 


3.3.  Prnyramn^ 

Copies  of  snomo.f,  snoao.daC,  decid.dat ,  conif.dat,  mix.dat, 
lapse. daC  and  SNOHO.RES  are  available  on  the  floppy  disc  Inserted  at 
the  back  of  this  volune.  A  fully  commented  hardcopy  of  the  programme 
code  Is  presented  In  the  following  pages. 


m 


j 
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C - - - - 

C  SMOMO.F,  Jim  1991 

BY 

KATBntin  K.  SAMBUS 

C . . . . - - - - 

c 

C  IBIS  IS  A  DISnUBlRBD  mSICAUT-BASED  SNOMELT  HDDEL.  SMONO.F  IS 
C  DBSIGHED  TO  OFBIATE  11  COmUKTICM  WITH  A  6IS  OECISIOH  AMD  DATA 
C  MAMIFinLATIOH  STHUCTURB  WHICH  DIVIDES  THE  CATCHKEMT  THAT  IS  TO  BE 
C  MODELLED  IRTO  BCMOCTmOUS  SUBDIVISIONS,  CALLED  CELLS.  SMOMO.F 
C  CALCULATES  THE  EMEBGT-BUDGET  OF  THE  SHOWPACE  AT  THE  HIO-FOIHT  OF 
C  ORE  CELL  AMD  CALCUUIES  THE  SESULTAMT  CHAN^  IM  SMOWPACX  nOFERTIES 
C  AMD  meltwater  RUNOFF.  USING  THE  RESULTS  FOR  EACH  CELL  AMD  THE  GIS 
C  STRUCTURE  A  MAF  OF  THE  SMOHCOVER  DISTRIBUTION  AMD  DEPTH  DISTRIBUTICH 
C  CAM  BE  OBTAIMED.  SMOMO.F  IS  DEVELOPED  FROM  THE  WES  MODEL  TSTM  AMD 
C  THE  SUBHUEL  VB6IE  HAS  IMOQRPORATED  IRTO  SMOMO.F  IM  MARCH  1989  WITH 
C  THE  HELP  OF  DR.  RAMDY  SCOGGINS,  WES. 

C 

C  SMOMO.F  IS  DIVIDED  INTO  2  MAJOR  PARTS; 

C  (1)  THE  MAIM  PROGRAM.  THIS  CALCULATES  MELT,  MANIPULATES  THE 
C  SROWPACK  DEPTH  AMD  PHYSICAL  CHARACTERISTICS  AMD  HANDLES  THE 
C  MODEL  INPUT  AMD  OUTPUT. 

C  (2)  THE  SUBROOTIME  TSTM.  THIS  CALCUUTES  THE  COMPONENTS  OF  THE 
C  EMERGY-BUDGBT  OF  THE  SHOWPACE. 

C  SMOMO.F  REQUIRES  THE  DATA  FILES  SMOMO.OAT,  LAPSE.OAT,  OECID.OAT. 

C  MIX. OAT  AMD  CONIF.DAT  IM  ORDER  TO  OPERATE.  THE  OUTPUT  FILE  IS 
C  SRQMO.RES.  A  FULL  LISTING  OF  ALL  THE  VARIABLES  USED  IS  INCLUDED 
C  AT  THE  END  OF  THE  PROGRAM  CODE. 

C  THERE  ARE  3  CODE  OPTIONS  FOR  SHOMO.F; 

C  OPTION  1  IS  THE  DEFAULT  OPTION.  OPTION  1  IS  UTILISED  UNLESS  OTHERWISE  STATED. 
C  OPTION  2  GIVES  AN  OUTPUT  OP  HOURLY  ENERGY-BUDGET  VALUES. 

C  OPTION  3  INCLUDES  A  MASS  CORSEXVATION  ROUTINE. 

C 

C - 

c 

C  IRITIALISATION,  DECLARATION,  DATA  STATEMENTS. 

C 

C0M«]R/MATRX1/XXX(  30 , 10  ) ,  TTY  (  30 , 10  ) 

OOHNNI /MATRX2 /DEPMC 1  ( 2 , 2 ) ,  DEIME2  (  3 , 2  ) 

C0H«]R/MATRX3/TBIM(1  ( 1 , 6) ,  TBIMa(2 , 6} 

O0»»EM(/VECT0R/VMENBN(  1 , 6) ,  VOLOSNd,  8) 

COMCN/TSTHl  /  PRESS ,  NCLOUD ,  ZA ,  SLOPE  1 ,  SURFC 1 ,  DAY ,  LAT 

C0m0R/TSlM2/H0MATL ,  RIPTS 

COM«MI/VEOI/JVEO,SI6F,8TATB,EPP,ALBEDO,BPOL 

COPKH/RAIR^AIRTP 

C0M«XI/SPItQ/SFItQ(6) 

SEAL  SRPFTH,8NVX,CRITDP,BCAPSR,DERSM, 

A  OERH.BCAPH,LBBATF,DEFTB,FTERM,RETRAO,AVSTP, 

A  AVBTT?,PFTN,DAYl,LAT,TGnR,TSOL,TABSOR,TATERM,TBTERM, 

A  TDTXBM,VSOIL(l,6),MBl,nDmi,HMB,l«8N,BOnT,ST?.E8NT?, 

A  SPBTI,SPRIH,SHDP1,TMAX,TM1M 
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nmoa  I.J.a,Kl,K2.Ka,IC10.K30.KOirntl..ILAFSE.IVEG,ISIIOH, 

A  OlAia.CZLL.TBAit 

c 

C  LATEMT  BEAT  OP  PUSIOB  IS  0.334 

C  SPECmC  BEAT  OP  ICE  IS  OIVEB  IB  •  2.10E+3 

C  SPECIFIC  BEAT  OP  HATER  IS  6IVEB  IB  -  4.1BE43  Ag'^'K' 

C 

DATA  OEIIH,aCAIW,BCAI>SB,LBEATP/1000.0.4.21,2.09,0.334/ 

DATA  SFBTI,SFBTH/0.0021,4.18/ 

DATA  TOTTIM,TPIQ.TPIIltT/l. 1.0, 60.0/ 

HIH-IO 

BCOHIP-ll 

BLAPSE-U 

RDECn^lS 

IHIX>14 

aOUT>9 

OPEHCUBIT-RIR.PILE-'SIIQHO.DAT*  ,STATUS-’OLO' ) 

REHin)  RIB 

OPEB(UBIT^ICOBIP.PILE— COBIP.DAT' .STATUS- 'OLD' ) 

REHIBD  BCOBIP 

OFEH(UBIT-BUPSB,PItE-‘UPSB.DAT'  .STATUS- 'OLD' ) 

REHIBD  BLAPSE 

OPEBCUHIT-RDECID.PILE-'DECID.DAT* .STATUS-’OLD' ) 

REHIBD  BDECID 

OPEB(OBITHHIX,PILE-‘MIX.OAT'  ,STATUS-’(EJ>* ) 

REHIBD  MX 

OPEB  ( UBIT-NOOT ,  FILE- '  SHOO .  RES  ’ ) 

REHIBD  BOUT 
C 

READ(10,*)CELL,TEAR 

C 

C  VECTORS  VBEHSB,  VQUSR  ARD  VSOIL  COHTAIR  TEE  FBTSICAL 
C  ISOPERTIES  OP  REH  8R0H,  CKJ)  SBOH  ARD  SOIL  RESPECTIVELY. 

C  VECTORd.D-  0.0 

C  VECT0R(1,2)-  SBOH/SOIL  TBEBMAL  DIFFUSIVITY,  ca^mln'^ 

C  VECTORd.S)-  SBOH/SOIL  BEAT  CORDOCTIVIIY, 

C  VEC10Rd,4)-  SBOH/SOIL  MSSIVITy,  DECIMAL 
C  VECTORd.S)-  SBOH/SOIL  ALBEDO,  DECIMAL 
C  VECTQRd.6)-  SBOH/SOIL  DEBSITT, 

C 

READdO,*>((VBEHBHd,J),J-l,6>,I-l,l) 

READ (10 , * ) ( ( VOLOSB (I , J } , J-1 . 6 ) , I-l . 1 ) 
READdO,*)((VSOIL(I,J),J-l,6),I-l,l) 

READ (10 ,*) PRESS , RCLOOD . ZA 
BEA0d0,*)8U»El,SURFCl,DATl,UT 
READ(10,*)ILAPSB 
READ(10,*)1VB6 
READ(10,*)Sia)Pl 
C 

C  IVIB  DEIBIHIHBS  MBICB  VBORATICH  DATA  PILE  IS  USED; 

C  IVEO-l,  OOnP.DAT  (OORIPBROOS  DAM) 

C  miM,  DECID.DAT  (DECIDUOUS  DATA) 


C  mo-3,  HU. DAT  (MUEO  DATA) 

C  XVEOO,  NO  VB3TATI0N  DATA  FIU  XS  USED,  THE  VEGETATION  IS 
C  MODELLED  AS  FASTUHE  (OPEN)  AND  CLEABCUT. 

C  ILAFSE  DETERMINES  HBICB  ELEVATION  DATA  FILE  IS  USED  TO 
C  HXIFY  THE  AIR  TEMPERATURES  WITH  LAPSE  RATE.  THE  FRESERT 
C  USE  OF  ILAPSE  REFERS  SOLELY  TO  THE  APPLICATION  OF  SIKRO.F  TO 
C  THE  H3  CATCBOn. 

C  ILAPSE-0.  ELEVATION  RAND  1300-2000FT.  THIS  IS  THE  ELEVATION 
C  BARD  (ff  THE  WTEOROLOGICAL  STATION  THAT  PROVIDES  THE  INPUT 
C  DATA,  THEREFORE  HO  CHARGES  ARE  MADE  TO  THE  AIR  TEMPERATURES. 

C  ILAPSE-1,  ELEVATION  BAND  1000-1300FT.  THIS  INCREASES  THE 
C  MAX.  AIR  TEMP.  BY  2’F  AND  THE  MIN.  BY  1*F. 

C 

IF(IVEG.BQ.1)THEH 

READ  ( 1 1  ,*)  SIGF ,  STATE ,  EPF ,  ALBEDO ,  HFOL 
ELSE 

IF(IVEG.EQ.2)THER 

READ  ( 1 3  .*)  Sm ,  STATE ,  EPF ,  ALBEDO ,  BFOL 
ELSE 

IF<IVBS.EQ.3)THER 

READ( ,* )SIGF , STATE , EPF , ALBEDO , HPOL 
ERD  IF 
ERD  IF 
END  IF 

IF ( ILAFSE. 6T.0)THEH 
READ(12,*)TMAX,IMIN 
END  IF 

RXAD(10,*)CRITDP 
READ( 10 , * )SOILOP, SOILTP, SSOLTP 
READ(10,*)N 
C 

C  N-  NUMBER  OF  DAYS  TO  BE  HXELLED  IE.  NO.  OF  DAYS  IN  DATA  FILE 
C  MINUS  ORE. 

C  DAILY  METEOROLOGICAL  DATA. 

C  XXX(1,8)-  JULIAN  DATE 
C  XXX(1,1)-  (nSERVATION  TINE,  24HR  CLOCE 
C  YYYd.l)-  MIRIMUH  AIR  TEMPERATIDIE ,  *C 
C  YTT(1,9)-  MAXIMUM  AIR  TEMPERATURE,  *C 
C  yYT(l,2)-  RELATm  HUMIDITY,  2 
C  YYY(1,3)-  CLOU^  COVER,  0-1 
C  YyY(l,e)-  HIND  SPEED,  ■/• 

C  YTY(1,7)-  PRECIPITATION,  an  HATER 
C 

RXAD(10,*)XXX(l,8),XXX(l,l),TT¥(l,l).YYY(l,0),YyY(l,2), 
A  TYY(1,3),TTY(1,8),TYY(1,7) 

READ(10,*)XXX(2,8),XXX(2,l),TYY(2,l),YYT(2,8),YyY(2,2), 
A  TYT(2,3),TTT(2,6),TTT(2,7) 

C 

C  SIMPLE  MODIFICATION  OF  AIR  TEMPERATURES  FOR  LAPSE  RATE 
IFdLAPSB.RQ.OJTBBR 
0010  S 
ELSE 
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C  COKVERI  ‘C  TO  ‘P 

m(i,i)-((m(i,i)*9.o)/s,o)+32.o 
m{l,9)-((YYY<l,8)«9.0)/S.0)+32.0 
YYY(2.1)-<(m(2,l)*9.0)/S.0)+32.0 
m(2,9)-(<YYY<2.9)*9.0)/5.0)+32.0 
IFdUFSE.BQ.DTHEH 
YYY(l.l)-m(l,l)+7MIK 
m(i,9)-yn(i.9)'fn«x 
YYY<2 , 1  >-TW(2 . 1  )+TMIII 
YYY(2.9)-m(2.9)+TM«C 
ELSE 

yifY(  1 , 1  )-TTT(  1 , 1 ) -IMIll 
YYy<i.9)-ma,9)-TM«c 
yYy<2,  i)“TW<2,  D-imH 
YYY(2,9)-YYT<2,9)-TMMC 
EED  IF 

C  CONVERT  *F  TO  ‘C 

m(l,l)-(<m(l,l)-32.0)/9,0)*5.0 
YYY(l,9)«<(TYY<l,9)-32.0)/9.0)*5.0 
yYY(2,I)-(<TTT(2,l)-32.0)/9.0)*5.0 
TYY(2.9)-<(YTY(2,9)-32.0)/9.0)»5.0 
END  IF 
C 

C  THE  SROHDEPTH  AND  THE  PHYSICAL  PBOFERTIES  OF  THE 
C  SIKWPACK  ARE  HELD  IH  THE  MATRICES  DERKl,  DEFMiU 
C  (HOLDING  THE  SROHDEPTH  AHD  OORRESPORDIIiG  TEMPERATURE) 

C  AND  TBKMd,  THXMXE  (BOLDING  THE  TBICXNES5  OF  THE  LAYERS 
C  OF  THE  PACK  AHD  THE  CORRESPONDING  PHYSICAL  PROPERTIES, 

C  IE.  VNEHSN  OR  VDLDSN).  THE  NUMBERS  1  AND  2  REFER  TO 
C  EITHER  A  1*  OR  2-LAYERED  FACE,  FOR  A  1-LAYERED  PACK 
C  NOMATL-l  AHD  NIPTS-2,  FOR  A  2-LAYERED  PACK  N0MATL>2, 

C  NIPTS-3. 

C  SET  2-LAYER  SNOH  THICKNESS  MATRIX  TO  ZERO 
C  SET  2-LAYER  DEPTH  MATRIX  TO  ZERO 
0 

S  DO  20  1-1,2 

DO  10  J-1,6 
THKMX2(I,J)-0.0 
10  CORTIRUE 

20  CONTINUE 

DO  AO  1-1,3 
DO  30  J-1,2 
DEIMX2(I,J)-a.O 
30  CONTINUE 

AO  CORTIRUE 

C 

C  SET  CONTROLS  TO  ZERO 
C  HRITE  HEADER  FOR  RESULTS  FILE 
C 

KORTRL-O 

SNVOL-O.O 

OEPTIi-O.O 


Kl-0 

K2-0 

K2S-0 

MRITE(9,*)'CEU.  MO.:', CELL 
WUTEtS, •} ’yBM;  ■  .yZAX 

HRITE(a,*)'JULIiUI  SMOHDEPTB  SMCMMELT 

&  DEHSITY  HAIM-CR-SNOH* 

NRITE(9,*)’  DATE  (CH)  (CM  SHOW)  ((M  HATER) 

A  (6/M**3)  (l-TES,  O-HO)' 

C 

C  DAILY  CALCULATIOM  LOOP 
C 

DO  3001  L-1,H 
Kl-Kl+1 
K30^ 

EELV1M-YYY(  1 , 1  )+273 .  IS 
PPTM-YYY(1,7) 

AIRTP-YYY<1,1) 

DDATE-XaCd.O) 

TIIC2-X)aC(2,l} 

C 

C  LOGIC  STRIICTUIIE  TO  IMITIALISE  SROHDEPTH 
C 

IP(IC3.EQ.1}IHER 

IBM 

GOTO  37 

Kl.gg 

IF ( PPTM . 6T . 0 . 0 . AMD . AIRTP . LT . 0 . 0 )TBEH 
SMPPTH-( PFTH/VHEHSM ( 1 , 6 ) ) / 10 
END  IF 
EHD  IF 

IF(K1.BQ.1)TBEM 

K2-0 

SMPPTII««in)Pl 
EHD  IF 

DEPTB-SRFPTM+DEPTB 
IF(DEPTB .GE .CRITDP)THEH 
IF(DEPTa  .GT  .CRITDF)THEI( 

ROHATL-E 

RIPTS-3 

DEIMX2(1,1}>0.0 

DETN(2(1,2)-AIRTP 

OEIMX2(2,l)-S.O 

DERa2(2,2)>AIRTP 

DEIMa(3,l)-SRPPTR 

DEma(3,2)-saoLTp 

XSRTP»((0EINC2(1,2)'*’0EPM(2(3,2>>/2)4’273.1S 
DO  S3  J-1,6 

TBM(2(1,J)^VRBMSR(1,J) 
nEMa(2,J)^V0LDSR(l.J) 

OORTIRUE 
THEMQd.D'S.O 


SHOHFALL 

(CM) 


S3 
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TI!I9«2(2 , 1  )>8IIFPTII-S .  0 

ELSE 

KMATL-l 

nFTS>2 

DER«1(1,1)>0.0 
DEEMK1(1,2)>AISTP 
0E1MC1(2,1)>0EPTH 
OEIMCl  (  2 , 2  )«SSaLTP 

ESIITP><  (DBFMCl  ( 1 , 2  )>OEfMEl  (2 . 2  )  ) /2  H273 . 13 
DO  34  J-1.8 

TBXMEl  ( 1 ,  J  )-VHEMSll(  1 ,  J  ) 

34  CORTmUE 

THia«l(l,l)4EPTB 
EHD  IF 

ISSa^ISMOWfl 

K2S>0 

K3>1 

E2«l 

eon  37 

ELSE 
MMATL-l 
EIKS-2 
DO  36  >1,6 
TBKMCK 1 ,  J>VSOtL(  1 ,  J) 

36  OORTIIIUE 
TBENCKl.D-SOILDP 
DEINCKl.D-O.O 
DEMC1(1,2)«AIKTF 
DE»«1(2,1)-S0ILDF 
DEmXl(2,2)-S0ILTF 
E3^ 

K2M) 

EHD  IF 
C 

C  CALL  SUBROOTIIIE  CSFUQ,  TO  CALCULATE  SFRQ 
C 

37  CALL  CSFRQd.S) 

C 

C  SET  DEPTH  MATRICES  AMD  VARIABLES  n  EERO 
C 

IFdKMATL.EQ.  DTBEH 
DO  SB  1-1,3 
DO  36  >1,2 
DEIM(2(1,J)-0.0 

38  canuniB 

39  OORIHDE 
ELSE 

DO  61  1-1,2 
DO  60  >1,2 

rapi«i(i,j)-o.o 

60  OOMTimil 

61  OCRlROE 


EUD  IF 
K1(M) 

AVSTP-0.0 
STP-0.0 
AVS(X.-0.0 
AWaR-O.O 
AVABSK-O.O 
AVATDt-0.0 
AVHTER<^.0 
AVDTEIt>0 . 0 
AVBTTF-0.0 
BOTTF-0.0 
C 

C  CALL  SUBBODTIIR  TSIM,  TO  CALCULATE  THE  SBORTHAVE,  LOKGHAVE,  LATEHT 
C  HEAT  AHD  SEHSIBLE  HEAT  EXCHAMSES  AT  THE  SURFACE  OF  THE  SHONFACK  AMD 
C  HEAT  TRARSFORT  THROUGH  THE  FACE.  IR  SRCMO  TSTH  CAR  MODEL  THESE 
C  EXCBARGES  FCR  EITHER  SROH  OR  SOIL. 

C 

CALL  TSTHd,  'T* ,  'R'  ,K10,AVSTF,AVSOL,AVraR,AVABSR,AVATER, 

A  AVBTER,AVDTER,K2,AVBTTP,TraR,TSOL,TABSOR, 

A  TATERM.THTERM.TDTERM) 

C 

C  FROGRAH  CORTIRUES  IP  SROH  HAS  MODELLED  BY  TSTM,  JUMPS  TO  THE  HEXT  DAY 
C  IF  SOIL  HAS  MODELLED 
C 

IF(E3.Eq. OTHER 
GOTO  9000 
ELSE 
K3-1 
ERD  IF 
C 

C  CALL  SilBROUTIRE  CMELT,  TO  CALCULATE  THE  MELT  RATE.  THE  ERERGY 
C  IRTRODUCED  TO  THE  PACK  BY  RAIR  (PTERH)  IS  HOT  IHCLUDED  AT  THIS  STAGE 
C  MELT  RATE  IS  CALCULATED  IH  CM  OF  SROH  ARD  Mt  OF  HATER  EQUIVALERT. 

C 

FTEMM.O 

CALL  aCLT(I,H,RETRAD,FTEIM,TSOL,TCTR,TABSOR, 

A  TATEiM,THTERH,TDTERM,PFTH,DERSH,i«ME,MISR,6TEm,DERH,DEFTH, 

A  AV8TF,AVBTTF,X8HTF,SFIITI,K30) 

IRAIH-0 

C 

C  CALL  SUBRODTIHE  RAIMIF.  TO  CALCULATE  THE  ERERGY  IRTRODUCED  TO  THE 
C  MELTIHG  PACK  BY  RAIR,  IF  RAIR  OCCURED.  THE  MELT  RATE  IS  RECALCULATED 
C 

IF  ( FFTR .  OT .  0 . 0 .  ARD .  ADtTF .  GT .  0 . 0 )  THER 
IFdfME.GT.O.  OTHER 
IRAlR-1 

CALL  RAIRMF(I,R,FFTR,FTERM,0ERH,X3O 

CALL  OCLT(I,R,RBTRAD,PTRRH,T80L.TQBR,TABSOR, 

A  TATIRM,nmgH,TDTBRM,FPTR,OERSR.M»B.lflSR,aTEiM,DBRH, 

A  OEFn.AFSTF.AVBTTF.XSRTF.BFBTI.KSO) 

BHD  IF 


no  IF 


c 

C  CAICUIATICM  Ot  SROUPALL  IF  MOT  FIBST  ITESATIOH  OK  SOIL  HOOELLED 
C  IM  TBE  FSEVICUS  ITBBATIOM. 

C 

IFIATSTF.GT. 0. 0)I8EM 

AV8TF-0.0 

ESD  IF 

IF(K2.BQ.0)THER 

IF(l>PIM.6r.0.0.AllD.AIRTP.LT.0.0)TBEH 

SMPPTM-(PPTH/VIIEHSH(1,6))/10 

X23>0 

ELSE 

SMPPTH-O.O 
E2»«2S+1 
EHD  IF 
EMD  IF 
C 

C  LOGIC  STROCTIffiE  TO  ALLOW  FCK  THE  EFFECTS  (ff  SMOWFALL 
C 

tF(Kl.BQ.l}TBEM 

OEiMC2<l,2)>AVSTP 

DEIMU(3,2)«4kVaTTP 

OEiW2<2,2)-(<DEmC2(2,  l)/OEPTH)*(AVBTIP-AVSTP)  )+AVSTP 
GOTO  3000 
EHD  IF 

IFCSHPPTM.GT . 0 . 0 )THEH 
IFCNOMATL.BQ.DTHEH 
IF ( TBENC 1 ( 1 . 2 ) . EQ . VMEWSM ( 1 , 6 ) ) THEM 
IF(ISinW.OT.O)TBEH 
TBEMCl  ( 1 . 1  )-TBIM(l  ( L .  1  HSHPFm 
OEFTB>OEFTB+SirPFTM 
DERai(2,l)-DEFia 
DEIM(l(l,2)iAVSTP 
DEmEl(2.2)-AVBTTF 
GOTO  4000 
ELSE 

ME1-(TBIKC1(  1, 1)*10)*TMEHSH(1,6) 

DEFTE1-(HE1/VQU>SH(1,S)>/10 

HaMATL-2 

HIPTS-3 

DO  70  j-i,e 

TBEMQC 1 ,  J)^VHBWBH(  1,  J) 

TBia«2(2,J)^VDL08H(l,J) 

70  OOMTIRUE 

TBEN(2(1,1)-8MPFTH 

THENU(2.I)-DEPTH1 

DEPTH-5HPPTH4DEPTH1 

DEII«2(1,2)^V8TP 

DBn«2(3,2)^VBTTP 

OBIHa(3,l)-DBPTB 

DEIHa(2 , 1  )-TID«2(  1 , 1) 
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Draa2(2 , 2)>(  (DKFI«2(2 , 1) /DEPTB)*(AVBnP-AVSTP )  )-fAVSTP 
OOTO  3000 
Dn>  IP 
ELSE 

l»MTL-2 
MIPTS-3 
DO  80  J-1,6 

TBItMC2  ( 2 ,  J ) -VDLOSM  ( I ,  J ) 

COMTIHUE 

THIHUd.D-SlIPPni 

TBKME2(2,1)4EPTH 

OEPTMEPTB-fSIlPPni 

DEFia2(l,2)-AVSTP 

0EmE2(3.2)*AVBTrP 

DEIH(2(3,l)-OEPIH 

DEIMC2(2 , 1  )*TH»«2(  1 , 1 ) 

DBR«2(2,2)-(  (0EPMX2(2.  l)/DEPTH)*(AVBnP-AVSTP)  )+AVSTP 
GOTO  3000 
END  IP 
ELSE 

IF(ISKM.(lt.O)TBEH 

TBEMX2( 1 , 1 )-THKHX2( 1 , 1 HSHPFTH 

DEPTB>DEPrsfsirppni 

DEFMU(3,1)-DEPTH 

DETN(2(2, 1)*TBI!MK2(1, 1) 

DEIMC2(2,2)-AVSrP 

DEmX2(3,2)>AVBTTP 

DEmC2(2,2)-((DEmC2(2, 1)/DEPTB)*(AVBTTP-AVSTP)  )+AVStP 

OOTO  3000 

ELSE 

HEl-(TH»gC2(l,l)*10)*ViaHSIf(1.6) 

DEPTai-(HEl/mDSII(l,6)  )/10 
TBia«2(2 . 1)>TBIH(2(2 , 1  )-H)EPTHl 
TflIM(2(l,l)-SIIPPTR 

DEPTB-THn«2  ( 2 , 1 ) '^TBEMU  ( 1 , 1  > 

I>EmC2(3.1>-DEPTB 
DEPM(2(2.  I  )-TBI9a2(  1 , 1) 

DEIMC2(1,2)<41VSTP 

DEPM(2(3,2)-AVBTTP 

DEIKt2(2 , 2)-(  (DER«2(2 . 1 )  /DEPTH)*(AVBTTP-AVSTP)  )+AVSTP 
GOTO  3000 
ERD  IP 
ERD  IP 
ELSE 

IP(R0HATL.EQ.1)TIBR 
IP ( TBEHCl ( 1 , 2 ) . IQ . VREH8R ( 1 . 6 » THER 
IPdSROH.Oi.DTHBR 
0EIMC1(1,2)>AVSTP 
DEM(1(2,2)-AVITTP 
GOTO  *000 
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0EPIB1-(HE1/«XI)SII(  1 . 6)  )  / 10 
OEPTB-OEFTHl 
DO  230  J-1,6 
TBMCl  ( 1 .  J  )-VQL08«  ( 1 .  J  ) 

230  OORTUniE 

TBKHKKl.D-OEPTH 
DER«1(1.2)-AVSTP 
0EPHC1(2,2)-AVBTTF 
0ER«1(2,1)>0BPTB 
GOTO  4000 
DID  IF 
ELSE 

0EEMEl<l,2)-AVStP 
OEmCl  (2 , 2  )-AVBTTP 
GOTO  4000 
EMD  IF 
ELSE 

IF(ISIiaH.6T.0)TBEII 

DER«2(1,2)-AVSTP 

DEFM(2(3,2)-AVBTTF 

OEINC2(2,2)-((DBINC2(2,l)/DEPTB)*(AVBTTP-AVSTP>>>AVSTP 
GOTO  3000 

njfft, 

NE1>(TBEMK2(  1 , 1  )*10  )«VIIEHSil(  1 , 6) 

DEPTH1>  (HEl /mOSX  1 , 6  )  )  / 10 
OEPTS-DEPTBUTHEMK2  ( 2 . 1 ) 

KHATL-l 
RIPTS-2 
DO  330  J-1,6 
TBKMXl  ( 1 ,  J  )-VDLDSI((  1 ,  J  ) 

330  oaiTiiniE 

TBENCKl.D-OEPTH 
DEIM(l(l,2)-AVStP 
DEI»«1<2,2)-AV8TTP 
DEn«l(2,l)-DEFTB 
GOTO  4000 
END  IF 
END  IF 
END  IF 
C 

C  LOGIC  STKOCTURE  TO  ALLOW  FOR  THE  EFFECTS  OF  SHOWELT 
C  ON  A  2-LATER  PACK. 

C 

3000  IF(DEPTB.LE.3.0)THEN 
GOTO  3002 
END  IF 


C  OPTION  3-MASS  C0R8BITATI0M  ROUTINE. 

C  IF  OPTION  3  RIQDIIID  INSERT  MASS  CONSEXTATION  ROOTIRB  TO  REPLACE 
C  ALL  CODE  HETHBEN  THE  ••••****LIHES . 

C 
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3001 


«00 


3500 


C 

C  LOGIC 
C  OH  A  : 

c 

AOOO 


4001 


iFHHGn.oB.DEFnima 

GOTO  3002 

ELSE 

1F(MRSH.LE.0.0)TBEH  • 

DEn«2(1.2)«AVSTP 

DEna2(3,2)-AVBTTP 

DEIW2(3,1)-DEFTB 

DEIMES  ( 2 , 1 )  •mOflU  ( 1 , 1 ) 

OEROE  (2 . 2  )-(  (DEIMK2  (2 , 1 )  /DEPTH)*  (AVBTTP-AVSTP)  l-fAVSTP 
TB»«2(2, 1)-DEPTB-TBD0(2(1, 1) 

GOTO  3300 
ELSE 

IF  (MISH .  GE .  TBIEMC2  ( 1 , 1 )  )THER 
■OHATL-l 
HIPIS-2 
DO  400  J-1,8 
IBEHCK 1 .  J)-VOLDSM(  1 ,  J) 

COHTIHUE 

DEPTB-OEPIB-MRSH 

nilMCl(l,l)>DEFTB 

DEIM(1(1.2)-AT8TP 

DEn«l(2,2)««VBTtP 

DEimKS,  1)>TBID«1(1, 1) 

ELSE 

TBEM(2  ( 1 , 1  )-Tn9«2  ( 1 , 1 ) -mSH 
DEFIH-TBEH(2  ( 1 . 1) -fTHEMES  ( 2 . 1 ) 

0EPMC2(1.2)-AVSTP 
DEIia2(3,2)-ATBTTP 
OEItSC2<3.1)>OEPTB 
DEPME2(2,  l)-THiafiC2(l,  1) 

DEPIV2(2,2)>(  (DER«2(2, 1)  /DEPTH)*(AVBTTP-AVSTF)  >+AVSTP 
EMO  IF 
GOTO  3000 
EHD  IF 
EHD  IF 

STROCTUHE  TO  ALLOH  FOR  THE  EFFECTS  OF  SHOWELT 
l-LATER  PACE. 


IF(DEPIB.LE.3.0)THEa 
GOTO  3002 
EHD  IF 

IF  (MISH .  a .  DEPTH)TBEH 

GOTO  3002 

ELSE 

IF(t«8H.LB.O.O)TBEH 
OEP»«l(1.2)^TSTP 
DEPHE1(2,2)-AVBTTP 
0EPMC1(2,  D-THEHCKl,  1) 
ELSE 

iffipnMnpiB-MiSH 

THMCKl.D-OEPTB 
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DERa(l(l,2)-AVSTP 
DEM(l(2,a)^VBTTP 
DEIMClCa,  D-TBXMCK 1. 1) 
m  IF 
EMD  IF 
C 

C  IF  TOP  LATER  OP  2-LATER  PACK  IS  <1CM  TBICIC  THEM  PACK  IS 
C  OSIVERTEU  ino  1-LATBS  PACK. 

C 

SOOO  IF(liaMATL.BQ.2.AirD.TBEKIC2(l,l).LT.1.0)THEH 

DEPTB-TBItt«2(  1,  D-fTHKMQCa,  1) 
imATL>I 

ifiprs-2 

DO  4SO  J*1,S 
TDMCK 1 ,  J)>VOLOSR(  1 ,  J) 

430  COHTIIIUE 

TBEMCKI.D-DBPTB 
DEPN(l(2,l)-OEPTB 
0EPI«1(1,2)-ATSTP 
DEPMIl(2,2)<iAVBTTP 
EHD  IF 
C 

C  MELTRATE  CAimOI  BE  BEGATIPE. 

C 

IF(tSSR.LT.O. OTHER 

MtSH-0.0 

EHD  IF 

IF(NIME.LI.0.0)TBER 

MBS-O.O 

EHD  IF 

330  FCmAT(6F6.2) 

333  FGRMT(2F6.2) 

C«****««***«**«*«*«**««****««*««****«**«******************«*********1 

C  BHD  OT  OPTIGH  3  CODE  IHSERTIOH. 

C 

C 

C  READ  IH  REXT  DAILT  DATA  LIRE,  SBIFTIHS  EXISTIH6  BOTTOM 
C  DATA  LIRE  TO  THE  TOP.  EXIT  PROGRAM  IF  LAST  DATA  LIRE 
C  REACHED. 

C 

9000  TYT(2,e)-TTT(2,6)/100.0 

TTT(2,2)-TTT<2.2)«100 
XXX(1,8)>XXX(2,S} 

XXX(1,1)-TII«2 

rTT(i,i)-rry(2,i) 

TTT(X.2)-TTT{2,2) 

m(l,9)-TTr(2,9) 

TTT(I,3)*m(a,3) 

TTY(l,8)-lfTT(a,6) 

tTT(1.7)-m(2.7) 

XXX(2,8)-0 
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m<2,i)-o 

m(2,8)-0 

m(a,2)-o 

m(2,3)-0 

m(2,«)-o 

m(2.7)-0 
WOOWHI-PlIiHI. 
irdDDQR.IQ.DTRBa 
OOTO  S002 
BHD  IF 

tinr<  1 , 1  )-TTY<  1 .  l)-273 .  IS 
YYY(l,fl)-TYY(1.9)-273.15 

BKAD(io.*)X]C]((2.a),xxx(2,i),yyy(2,i),m(2.8>,m(2.2), 
&  m(2,3).rfT<2.6),m(2.7) 
c 

C  CONVERT  AIR  TEKPERAIORE  FOR  LAPSE  RATE 
C 

IFdLAFSE.BQ. OTHER 

OOTO  90 

ELSE 

C  CONVERT  *C  TO  ‘F 

m(2,l)-<<trr(2.1)*9.O/5.0)+32.0 
m(2,9)>((Yn(2,9)*9.O/S.0)-t'32.0 
IFdLAFSE.EQ.DTBEH 
TTY(2, 1)-TTT(2,  D'dMIN 
TTT(2,9)-TYT(2.9HTHAX 
ELSE 

TTT(2 . 1 )-TTT(2 . 1 )*IMIN 
TTT(2,9)-TTT(2,9)-TMAX 
END  IF 

C  CONVERT  ‘F  TO  *0 

m<2,l>-((m(2.1)-32.0)/9.O*3.0 
m(2,9)«(<TTY(2.9)-32.0)/9.0)*3.0 
END  IF 
C 

C  REINITIALISATION  OF  CERTAIN  VARIABLES,  VOUTE  OUT  RESULTS. 

C 

90  IF(E2.EQ.0)TBER 

IF(SHPFTN.8T. 0 . 0 JTHEN 
ISN0M>ISN0W+1 
ELSE 
ISNOH-0 
END  IF 
END  IF 

AIRTP-TTTd.l) 

KELVm-TTYd ,  1  )-f273 . 13 
FFIN-TTT(1,7) 

D0AIB4XXd,S) 

BamLHBORm^i 


HRITB(0,333)DMTt-l,OBFTH,MISN,MWB.SN0H,DBNSN,  TRAIN 


i 
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333  r(aM*T(Fe.2.S](,P6.2,tX,r6.2.4X,P6.2.4X,rs.2,7X,F6.2,SX,I3> 

sool  oomimtt 

S002  MRXTE(9.*)'nOGSAM  TBBOaATIOII' 

STOP 

no 

c 

C . - - - - - - - 
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c 

C  SUBROOTin  CHELT 

C - 

C  THIS  CALCOLATIS  TBB  WLnATE  »  CM  QT  SKM  AK)  Mi  (7  KATER  EQUIVALEHT. 
C  THE  DSaSITY  Of  TIB  SKMPAOC  IS  CALCUUTED  IF  THE  PACK  IS  A  E-LATEBED 
C  PACK.  ICLT  IS  CALCOLAIED  OSIRS  THE  SHOW  EHEBGY-BUIXZT .  MOST  OF 
C  THE  OOWEHtS  OP  THE  SMOH  EHEaST-BUDGET  ARE  CALCULATED  BY  TSTM,  THE 
C  EXCEPTIOm  ARE  PTEBH  AID  DELQ6.  THE  EHERGT  OHITS  ARB 
C 

SUBROUTIHB  OCLTd.H.HETRAD.FTERM.TSOL.TOR.TABSOR, 

A  TATERM.THTERM.TDTEBH.PPTH.DEISR.MME.tCSR.STElM.DElIH, DEPTH, 

A  AVSTF,AVBTTP,ESITP.SPBTI,K30) 

C(Mai/HATIX2/DEII«l(2.2).DEfH(2(3.2) 

OMCII/HAIRXS/TBEMEKI  ,6)  ,TH»a2(2. 6) 

OOHCI/TSIME/KMATL  .  HIPTS 

REAL  LHEATF.HETRAD.MMB.MISI, DEPTH, DERSH.A.B.C.D, 

A  XSITPl, DELOS, DELIME, DELRSR,X8ITP 

DATA  LBEATF,BCAPH,BCAPSH/0. 334,4.21,2.09/ 

OTERM-O.O 

C 

C  CALCULATIOI  OT  THE  DEBSITT  OF  SIKM 
C 

IFdCMATL.EQ.DTHEH 

DEISIHTBEME1(1,6) 

ELSE 

A-TBEMa(l,l)/DEPIB 
B>THI(ME2(2,  D/DEPTH 
C-A*TBIME2(1,6) 

D-B*TB»«2(2,e) 

DEXSI-OO 
BID  IF 
C 

C  CALCULATIOI  OF  DELQS  II  THE!  (XIVERTED  TO  lUB~^dy'^ 

C 

IF(K30.BQ.1)TBEH 

E30-0 

OOTO  10 

ELSE 

XSRTPl^eiTP 

XSITF>( (AFSTP4AFBTTP)/2)>273 .13 
OELQB-((DE?IB*0.01)*(DEMSI*1000))*((SPBTI*XSITF)- 
A  (SPBTI*XSITP1)) 

BID  IF 
C 

C  CALCULATIOI  OF  THE  EHERST-BUDIXT  OF  THE  SIOMPACE  ABD  THE  MELT  RATE. 

C  MULTIPLY  MINE  BY  1000,  SO  THAT  RESULT  IS  IH  Wt  HATER. 

C  MULTIPLY  MMI  BY  100,  SO  THAT  RESULT  IS  IR  CM  SNOH. 

C 

10  mXAD-ISCL-CTSOL-TABSOD'f'TATERH-TIMR 

DEKRI  lEYRA&tlUHMIHUIERMmTIIMIfimM  DELOS 
UUM^dnLQH/  (LIBATF*DBIH)  )*1000 
mUWIDELOM/  (UBATP*  (10181*1000 )  )  >*100 
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c 

C  SUBROOTIHE  RAimP 
C . 

C  CALCULATES  THE  EMERST  IHPUT  TO  THE  PACK  BT  RAIN  FALLING  OH 
C  A  MELTING  PACK  AT  OOC.  THE  EQOATiaH  IS  TAKEN  FROM  MALE  A 
C  GRAY  (1881).  THE  UHITS  ARE  MJa'^dy'^. 

C 

SOBROCrriRE  RAIIMPd.H.PPTH.FTERH.CENW.XSO) 

COHCH/RAIN/AIRTP 

REAL  PTERM 

DATA  LHEATF.aCAPSN.HCAPW/0.334.2.09.4.21/ 

DATA  SPHTH/4.18/ 

PTERf^(DEHW*SPHTW*AIRTP*FFTH  )  / 1000 
PTERI^PTERM/1000 
K30-1 
RETURN 
END 
C 

C . . . 
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1 


1 

i 


C  SUBROUTIHE  CSFRQ 
C . - . 

C  CALCULATES  SFRQ.  THE  VERTICAL  GRID  SPACIIR3  IH  EACH 
C  COHPUTATIOHAL  LAYER. 

C 

SUBROUTIHE  CSFRQd.H) 

COMCW /MATRX3  /  THUff  1  ( 1 , 6  ) ,  TBEMC2  (  2 , 6  ) 

COMOH  /  VECTOR/VHEMSH  ( 1 , 6 ) .  VOL08R  ( 1 , 6 ) 
C0M«3H/TSTM2/HCMATL .  HIPTS 
COM5M/VEGI/  rVEG,  SICT ,  STATE .  EPF .  FOLA,  HFOL 
CCmON/SFRQ/SFRQ(  6 ) 

REAL  SFRQ1,SFR02 

DATA  TOTTIM.TFRQ,TPRHT/l, 1.0, 60.0/ 

C 

C  CALCULATION  OF  SFRQ  FOR  LAYER  1 
C 

IF(N0MATL.EQ.1)THEN 
SFRQl-CTHKMXK  1 , 1 ) ) /2 
ELSE 

SFRQl-(TBia«2{  1 , 1} ) /2 
EHD  IF 

100  IF(SFRQ1.LE.S.0)THEH 

GOTO  200 
ELSE 

IF(SFROl.U.10.0)TBEN 

GOTO  200 

ELSE 

SFRQl-SFRqi/2 
GOTO  100 
EHD  IF 
EHD  IF 

200  SFRQ(1)-SFBQ1 

C 

C  CALCULATION  OF  SFRQ  FOR  LAYER  2.  IF  NECESSARY 
C 

IF(HCHATL.EQ.2)THEH 
SFRQ2-THI(MC2(2 , 1 )  /2 
300  IF(SFRQ2.LE.S.0)THEN 

GOTO  400 
ELSE 

IF(SFRQ2.LE.10.0)THEH 

GOTO  400 

ELSE 

SFRQ2-SFRQ2/2 
GOTO  300 
EHD  IF 
END  IF 

400  SFRQ(2)-SFRQ2 
END  IF 


C 

C  CALCULATION  OF  TDCSTEP  EH8URIHG  NUMERICAL  STABILITY. 
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c 

C  SUBROUTINE  TSTM 

C - 

C  CALCULATES  THE  SBCitTHATE,  LOKUATE.  LATENT  BEAT  AMD  SENSIBLE  BEAT 
C  EXCBAN8ES  OVER  TBE  SMOHPACK  OR  SOIL  SURFACE,  AMD  TBE  BEAT  TRAMSKST 
C  TBROOGB  TBE  SOIL  OR  SHOW.  T8TH  BAS  BEEN  CONVERTED  INTO  A  SUBROUTINE 
C  FROM  BALICE  ET.AL.  (IRSla  A  b).  TBE  CONVERSION  OF  TSTM  AND  INCLUSION 
C  OF  VE6IE  HAS  AIDED  BY  DR.  RANDY  SCOGGINS  (WES).  VEGETATED  SURFACES 
C  ARE  MWELLED  USING  TBE  SUBMXEL  VEGIE.  TSTM  HAS  CRIGIMALLY  TBE  TERRAIN 
C  SURFACE  TD1PERATURE  MODEL  AND  ITS  FUNCTION  MAS  TO  PREDICT  SURFACE 
C  TIMPERATURE. 

C 

C 

C 

SUBROUTINE  TSTMd.ANl.AM.ElO.AVSTF.AVGBR.AVSOL.AVABSR, 

A  AVATER,AVBTER,AVDTER,K2,AVBTTF,IGBR.TSOL.TABSOR,TATERM, 

A  TBTERM.TDTERM) 

C 

C  INITIALISATION,  DECLARATION  AMD  DATA  STATEMENTS 
C 

INTEGZR  OUTCD 
ODCMSION  DATMK<10O,I0) 

DIMENSION  RBQC(6),MAX(10),DEm(*S0),FHH(30,10),BBB(30,10) 
OmENSIOR  TITLE(7) 

DDSNSION  MX(6),ATF(2),FEB(Z) 

DIMENSION  RR(6},IHTR(7) 

DIMENSION  TBE(6),ALFB(6),FK(8),STGR(7,4S0),5TABM(6> 

SEAL  KTaiF6,K7ENPA,LAT,ACL(8),BCL(S),M,ESQ,BOTTF,HET, 

A  ISaL,ICaR,IABSaR,IATERM,IBTERM,IDTERM,ISUSFG,IREFRA, 

A  SLOR,SURFAC,CLR(8) 

CBARACTER  BEADER*72,AR*1,AR1*1 
COMCN/MATRXl/EXXOO ,  10) ,  YTYOO,  10  ) 
COMON/MATKE2/DEF(V1(2,2),DEFMC2(3,2) 

CCMON/MATRXS/TBENEl  ( 1 , 8)  ,TBRM(2(2, 6) 

CCmON/TSTMl /PRESS ,  NCLOUD ,  ZA,  SLOFEl ,  SUSFCl ,  DAY ,  LAT 

C0»tai/TSTM2/N0MATL,RIPTS 

OOMCN /VEGI  /  IVE6 ,  S  1(7 ,  STATE ,  EFF ,  ALBEDO ,  BFOL 

CONOR/SFRO/SFRO(B) 

DATA  TDTTIM,TFRQ,TFRRT/1, 1.0, 80.0/ 

DATA  CLR/0. 04, 0.08, 0.17, 0.20, 0.22,0. 24, 0.24,0. 23/ 

DATA  OUTCD/0/ 

DATA  ACL/82. 2. 87. 1,32. 3, 39. 0,34. 7, 23. 8, 11. 2, 13. 4/ 

DATA  BCL/. 079,. 148,. 112,. 083,. 104,. 139, -.187,. 028/ 

DATA  8iaMA,FI,AC,BC/8.12E-U,3.141S93,17.268,33.88/ 

DATA  CC/0.2ei/ 

DATA  LAST, G,E8Q,CF/2, 980. 0,0. 18,0. 24/ 

C 

C  STATEMOIT  FUNCTIONS  FOR  USE  IN  VE(7TATIOR  SECTION 
C 

E(T)HIB*8. 108*EXF(AC*(T-273 . 13)/(T-BC) ) 

BSAT(T)-8. 10S*EXF(AC*(T-273 . 13)/(I-BC) ) 
Q(T)-0.022/(FRES8/E(T)-.378) 


z 
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qSAT(T)-0.622/(FBESS/mT(T)-.378) 

C 

c  convaaioH  of  variables  passed  psoh  maim  ncxsuM  imto  pcbm 
C  USED  n  TS1H.  1  CK  2-LATERED  SMOHPACE  USED. 

C 

IFLA62-0 

IRETRII-O 

IFIWT-O 

SUH-O.O 

IFdKMATL.EQ.DTBEH 

XXE(1,S)-DEIM(1(1,1> 

m(l,S)-DER«l(l,2] 

m(2,S)>DEFiai(2.1) 

m(2,3)>OEII«l(2,2) 

TBK(l)-<TaEMCl(l,l) 

ALPB(1)-TB»«1(1,2) 

FE(1)-TBEME1(1,3) 

EPSR-THKMEKl.A) 
aMLLA>(  1-TBIMEl  ( 1 ,  S ) ) 

ELSE 

XXX(1,S)-DER«2(1,1) 

m(l,5)-OEiMC2a,2) 

XXX(2,S)>DEm(2(2,l) 

m(2,S)-DEn«2(2,2) 

XXX(3,S)>OEm(2(3,l) 

nT(3,3)-DEIM(2(3,2) 

TBK(l)>TBngC2(l,l) 

TBK(2)-TBEME2(2,1) 

ALPB(1)-TBI9«2(1,2) 

ALFH(2)-T0KHK2(2,2) 

PK(1)-THEI«2(1,3) 

FE(2)-TBIEMC2(2,3) 

EPSII-THEME2(1.4) 

SMALLA^d-TBEHUd,  3) ) 

ERD  IF 

C - - - 

C  IRITIALIZE-VARIABLES-ARD-OCMSTAIITS 
98986  ASSIGH  9999A  TO  199993 
SO  TO  99893 
89994  CCMTIRUE 

C - 

C  IRPUT-OATA 

ASSIGN  89969  TO  199980 
GO  TO  99990 

C - - 

C  PRnrt-IRPUT-OATA 

99989  ASSIGN  99987  TO  199988 

IF(ANl.EQ.'T')aO  TO  99988 
99987  OOniRUE 

C - 

C  CALCULATE-TABLE-SLOPE-AND-IRTERCEPT 

99986  ASSIGN  99982  TO  199983 
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ao  TO  89963 

C - - - 

C  CAIjCULATE-SUIiF6CB-6ia>-LAyn-TBffEIlATUBE 
99962  CORTIRUE 

ASSIG9  99960  TO  199961 
00  TO  99961 

C - - 

99960  OONTimiE 

PPTII-<m(  1 , 7  )+YYY(2. 7 ) ) /2 
C 

C  CALCULATION  OF  DAILT  TOTALS  FOR  THE  EMEROY-BUDGET  VARIABLES 
C  REQUIRED  BY  THE  HAIR  FHOGRAH. 

C 

StMC2-0 

SOC3-0 

SUHCA-0 

SUHCS-0 

SUMC6-0 

SIMC7-0 

SUHCB-0 

SUC9>« 

DO  99  >1,K10 
SUMC2>SUNC2-H>AT(a(  1 , 2 ) 

SUMC3>«aMC3+OATI«  ( 1 , 3 ) 

SQMCA-SUMCA+OATMXC 1 , 4 ) 

SUC»4QMC»K)ATI«(I ,  S) 

SIMC6>«»Ce40AlNE(1 , 6) 

SUMC7-SUHC7«DATI«(  1 , 7 } 

SOCO-SOtCBfOAIMC  ( 1 , 6  ) 

SUMCO-SUMCB+DATNCd  ,9) 

99  COMTIHUE 
C 

C  CQRVERT  DAILY  EHEBOY  TOTALS  IH  lyhr'^  TO 
C 

TraR>SUHC3/24 

TSaL>SaC472* 

TABS0IMUMC3/24 

TATEBH4IM:6/2* 

THTBat-SUHC7/24 

TDTERH>SUMC6/24 

C 

C  CALCULATION  OF  HEAR  DAILY  EHERCY  VALUES  IH 
C  ARD  HEAR  SURFACE  ARD  BASAL  TZMFERATURES  (*C). 

C 

AVSTF-SaC2/K10 
AVOm-(SUHC3/E10  ) /24 
AV80L-(SaMC4/E10)/24 
AVABSie-i8IMCS/X10)/24 
AVATER^(S0HCe/E10 ) /24 
AVRTB»>(8»C7/E10  ) /24 
AVDTER-(8IMC8/E10  ) /24 
AVBTTM0MC9/K10 
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c 

C  RRUBII  TO  Moa  mOGBAM 

c 

miTB(*,*)‘mMkL  TCRHIMATiaR' 

C - 

C  STOP 

KETUSR 

C - 

C - - - 

99997  CORTnUE 
C  FOBHATS 

90  FOnMTC  BOTTOM  BOURltty  11IDBX-M3) 

92  POBMATC*  BOTTOM  BOUMIBIY  TEMPESATUBE-’ ,F6. 1,  ’  DEO  C*) 

93  FOBMAT(‘  BOTTOM  BOain»y  BEAT  FLUX-’, F6.1.‘H/H**2*) 

97  F0BMAT(F11.2,F11.1,F13.2,F13.1.FB.2) 

120  F0»MT(A8,F8.1) 

1*0  F0mAT(F3.1,F10.1,Fe.l.F12.1,F12.1,F13.1) 

US  FrBMAT<F8.1,F7.1) 

130  F0ilMAT(F12.1.I12,FU.2) 

160  FOIIMAr<F11.2.F13.2,F8.1.Fl0.1) 

170  F0IIMAT(F9.4,F12.1) 

180  FORMAT(I7,P16.1,Ftl.0,F12.0) 

C190  FaRMAT(*X,FS.2,F6.1,2X,t8,I9,lll,I9,IB,I8) 

190  FamAT(4X,FS.2,Fe.l,F8.2,r8.2,P8.2,F8.2,F6.2,Fe.2) 

260  FOnMAT(9X,F8.1,F11.2,F9.2,F9.2,FB.2.F10.2,F8.2) 

200  F0BHAT(F6.2.F8.2) 

210  F0BMAT(I6,F11.1,F12.1,FU.2,FU.2) 

220  FORMAT(lBl) 

230  FORMAT (A72) 

240  F0n<AT(AX,F7.2,3X,FS.0,SX,Fe.2,10X,F6.2,9X,F7.2) 

230  F0nNAT(6X,F6.1,13X,F*.l,12X,FS.l) 

310  FORMATCllX, 'TOTAL  GBATBGOT  EFFECTIVE  G80UBD  FOLIAGE', 
A  4X, 'SOLAR') 

320  FOmAT(UX,'RADIAllCE  TEMP’ ,10X, 'TIMP  TEMP' 

A  ,4X,’IIISOLATICM') 

330  FORMATISX,  'BX'  ,7X,  ’  (LAMGLETS)'  ,9X,  ’  (C)'  ,11X,  ’  (O' 

A  ,6X,’(C)  (LAHSLETS)’) 

3*0  F0mAT(SX,' - REFL-RREFL - REFL - RREFL' ,30(1B-) ) 

270  F0aMAT(3X,FS.2,F8.2,F8.2,2X,F8.2,F8.2,F8.2.F8.2,F8.2,F8.2) 

350  F0BMAT(S7X,’SEiBIBLE  LATERT') 

380  FOBMATCSX, 'BBS  SURFACE  GRATBOOT  SOLAR  SURFACE  ATMOS’, 
A  ’  IR  BEAT  BEAT') 

370  FCnMAT(llX,’TEMP  RAOIAHCE  IBSOLATIOH  ABSCRP  DaSSIOR’, 

A  ’  LOSS  FLUX’) 

380  F0BMAT(11X,'DB0  C ,23(1B-) , ' (LARQLEYS)' ,2*(1B-)) 

390  F0BMAT(3X,’TDC’,18X,'DEPTB,  TEMPERATURX’/*X, 'BR' ,21X, 'CM' , 

A  aX,’C'/2X,6S(lB-)) 

*00  F0BMAT(2B0  ,FS.2,*(3X,FS.l,’  ’,FS.2)) 

*10  F0BMAI(iaX.F3.1.1X,F5.2,3X,FS.l,lX,F5.2,3X,F5.l,lX,F5.2, 

A  3X,FS.1,1X,PS.2) 

GO  TO  199997 


C- 
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9999S  CCHTIIIUZ 

C  TO  mTULIZZ-ViOUBLES-iUn-CCIBTiUnS 
BB— 2.4E-* 
tUX(l)»2 
tUX(2)-2 
MM((3)>2 

HIUr(4)-2 

HM((3)«inPTS 

HAX(8)-2 

HAX(7)-2 

MAX{8)-0 

(MX(9)-0 

HM((10)>3 

IBOOO 

IBOP^ 

00  TO  199993 

C - - - 

99990  OOMTIXUZ 
C  TO  IBPUT-DATA 
C 

C  HffOT-ATMOSPanUC-SPBCIFlCATIOBS 

ASSIOB  99978  TO  199979 
GO  TO  99979 
C 

C  IBm-SURPACI-ORIEBTATIOB-SPBCmCATIOIJS 

99978  ASSIGR  99978  TO  199977 
00  TO  99977 
C 

C  nm-BEAT-rUM-CACaLATIOH-OQimiOLS 
99976  ASSIOB  99974  TO  199973 
00  TO  99973 
C 

C  IRFUT-IBITIAL-TIMFERAniEe-FItOFIU 
99974  ASSIOB  99972  TO  199973 
00  TO  99973 
C 

C  IHm-TOP-SURFACB-COBSTSBTS 
99972  ASSIOB  99970  TO  199971 
00  TO  99971 
C 

C  IBm-LATBK-SFBCiriCATiaBS 
99970  ASSIOB  99988  TO  199989 
00  TO  99969 

C 

C  IBFUT-BOTTOM-BOUBDART-DATA 

99968  ASSIOB  99986  TO  199967 
00  TO  99987 
C 

C  IBmT-VBOETATlOB-PAiUiCTBBS 

99988  ASSIOB  99798  TO  199799 
00  TO  99799 


C 
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8S7M  00  TO  199080 

C - 

89978  OOnimiE 

C  TO  maT-ATNOOnBRlC-SnCIPlCATIOB 
c 

c  lub  time  aib  temp  bb  clood  coveb  hird  speed,  ibsolaticm 

C  HO.  BB  DBG  C  S  (0-1)  M/3  CAL/CM**2-HIH 

c 

)Qa((2,l)>]aaC(2,l)'f2400*(XEX(2,a)-]QQ((1.8)) 

860  DO  8086S  J-1,2 

I»(X-I1IT(XXX(J,1)/100.0)*100 
Am-(XXX(J,1)  -  BEAL(QOQ())/60.0 
BEX  ( J .  1  )*«BAL  (  DQCX )  / 10  0 . 0>Am 
XEX(J,7)-XXX(J,1) 

XEX(J,2)-XXE(J.l) 

m(j.3)-m(j.i) 

XXE(J.4)-XXX(J,1) 

iaac(j.6)-m(j,i) 

TYY< J . 1)-TTT< J, l)+273 . 1 
YTr(J.9)-TTT(J,9)+273. IS 
YTT{J.2)-TTY(J.2)»0.01 
TTT( J. 6 )-TTY( J, 6 J*100 . 0 
TTT(J,*)-TTT<J.4)/697.6 
0808S  OOBTHIUE 

]aa((3.io)"]aa((2,i} 

]aa((2,iO)>i4.o 

XXX(1,10)«XXX(1,1) 

TrT<3,10)-TTT(2,l) 

TTT<2,10)-TYT(1,9) 

TTT<1,10)-TTT<1,1) 

C 

840  PACTH-(1000.0/PBESS)**0.286 
00  TO  188879 

C - 

99977  CGRTIiraE 

C  TO  IRPUT-SUBPACE-OUEBTATICXI-SPECIFICATIOliS 
C 

C  LIKE  SFC  SLOPE  SFC  AZSfJTB  DAT  LATITUDE 

C  BO.  DEB-BQRIE-O  DEO  S>0  JULIAR  DEG 

C 

SLOPE-SLOPEl 
SUBPAOSUBPCl 
SL(»E<4LOPE*PI  / 180 . 0 
SUBFAOSDBFAC*?!/ 180 . 

00  TO  189977 

C - 

9987S  OOniBUE 

C  TO  IBPUT-BEAT-PLOM-CACUUTiaR-OOBTROLS 
C 

C  LIRE  BO.  09  BO.  OP  24  BBS  TOC  STEP  PBIBT  PBEQ 

C  BO.  LATEX8  BBPITIOBS  MIR  MIR 

C  1-8  2-3 


c 


ioTTiM>]aac(2.i)  -  m(i.i) 

00  TO  100873 

C - 

80873  comma 

C  TO  IHm-inTIAL-TBVESAnU-FBOriLZ 

c 

c  Lin  K).  or 

c  K.  noriLE  niRTs 

c 
c 

C  Lin  OEPTH  TEMP 

C  HO.  CM  DEO  C 

c 

IUX(S)-HIPI8 
DO  00864  J3>1,HIPTS 
«Y(  J3 ,  S)-m(  J3 , 3)+273 . 13 
00064  OOHTIHUE 

00  TO  180073 

C . . 

80871  CQHTIHUE 

C  TO  IHm'TOP-SUKPitCE-CaMSTAirrS 
C  Lin  EmSSmTT  ABSOMTIVITy  HOISTOSE 

C  HO.  COHTERT  DEY  Wt. 

C 

FMrrA-SIGHA*EFSH 

c 

00  TO  108871 

C - 

88867  OQHTimiE 

C  TO  THFOT-BOTTOM-BOOHDARY-DATA 
C 

C  IF  LFLUn-0,  TnU  IS  HO  nAT  FLOX  TBROUGB  BOTTCM 

C  IF  IFLOXY  LT  0,  TBEU  IS  HO  AUtSFACE  BEHEATB  BOTTCM 

C  IF  IFLOXY  01  0,  THEn  IS  AHISFACE  BEHEATB  BOTTOM 

C 

LFLOXY— 1 

IF(.BOT.(LFLOXY.EQ.O))  00  TO  80062 
DI»IO-yYY(HlPTS,5) 

TB-0. 

00  TO  80863 

00862  IF(. HOT. (LFLOXY. LT.O)}  00  TO  08861 
C 

DFBMl— 0.3 
TB-FK(HaMTL) 

BEF-O.O 

BE-0. 

BEP-O. 

TIM.O 

FACTD-0. 

FACTE-O. 

n-0. 
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DnMi-DnBu/«87.e 
00  W  99S63 

98901  oGaniiias 

88803  00  to  X889e7 

C - 

98968  OOMTimB 

c  TO  mvT-um-snciFicATion 
c 

C  LIRE  TBICniESS  VERT.  GRID  TRERMML  OIFF  BEAT  OCMD 

C  ao.  CM  SPACE-CM  CM**2/HIR  CAL/MIR-CM-X 

C 

00  88900  JA'l.ROHAn. 

RaOC(  JA  )4K(  ) /ALFB(  J4 ) 

88900  OOMTIROE 

C  CORTBOL  SNITCHES  SPECIFIED 
C 

RIT*1 
IEFSHT-1 
00  TO  198900 

C - 

08798  CCRTIROB 

C  TO  IRPOT-VEIKTAnaR-PAXAWTERS 

C  IVE6  DETEBMIRES  MRETBER  VEGETATIOR  PARAMETERS  ARE  USED 
C  OR  HOT.  BEICBT  OP  THE  lOTBaROUGICAL  IRSTROteHTS  FROM 
C  NBICB  IRE  DAILY  DATA  ARE  TAXER  ARE  ASSUCD  TO  BE  lOOCM 
C  ABOVE  THE  FOLIA0E  HEIOBT,  SEE  TEXT  FOR  FULLER  EXPLARATIOM. 

C 

IF(IVBG.EQ.O)  00  TO  1120 

rOLA-l-ALBEDO 

IF(HF0L.a.2A)TBEH 

EA^FtX.AlOO.O 

ERD  IF 

IFCSIQF.LE. 0.0)00  TO  199799 
TF-TTT(1,1) 

IVBO-1 

EPl-EPPAEPSR-EPfEPSR 
ZO-0 . 131*HFaL**0 . 887 
CHOHCSQ/ (AL0G(ZA/Z0)**3) 

2DSP^ .  701*BFGL*^ .  878 
CaB-KSQ/(ALQO( (ZA-ZDSP) /Z0)**2) 

CaO>(l . -SiaF)*CHOASIQF«CBB 
DELIMIT. 

0AF-Q8AT(TF) 

1120  00  TO  180788 

C - 

98880  CCRTIHUE 
C  TO  PRIRT-IRPUT-DATA 
C  MUTEC^.ISO) 

MR1TE(*,220) 

HRITEI*,*)'  • 

MaTE(*,230)HEADER 
HUTRC*,*)'  ' 
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MutB(*,*)‘  Aneemuc-sPBcmcATioiB' 

WUTB(*.*)'  • 

Aim  FS£Sa  CLOOD  TYPE  SHELTOt' 

MUTBC*,*)'  W  IBKE  HEIOn-CM' 

MUTB(*,  lS0)mS8,KUIU0.ZA 
l«ltE(».*)'  ' 

MUTE(*.*)'  TDC  AIR  TEMP  RR  CLOUD  COVER  WIRD  SPEED' 

A  SOUR  IRRAO* 

VBITE(*,*)'  BRS  DBG  C  X  (0-1)  M/S  W/H**2’ 

(tMX-2 

l«m(*,U0)(lQQC(J,l),m(J,l)-273  1S. 

A  m(J.2)*100.0.¥YY<J.3),YYT<J,6)*0.01,m(J.*)*687.6.J-l.(*1M() 

l«ITB(*,*)‘  SORFACE-CSIBIITATXOR-SPECIPICATIOMS’ 
l«ITE(*,*)'  ■ 

MUTB(*.*)'  SPC  SLOPE  SFC  AZIMUTB  DAT  UTITUDE' 

t«ITB(*,*)‘  OEG-OaillZ-O  DEO  S-0  JULIAH  DEG’ 

MUTE(*,  160  )SLOPE*180/PI  ,SURPAC*180 . 0/PI ,  DAT,  UT 
MRITE(*.220) 

MIITE(*,*)'  HBAT-FLOH-CACCUTIOM-COinROLS' 

WUTEC*.*)'  • 

t«ITE<*,*)'  RO.  OF  NO.  (V  2«  BRS  TIME  STEP  PRINT  FREQ' 
MUTEC*,*)'  UTERS  REPEIITIORS  KIR  MIR’ 

WITEt*.*)’  1-6’ 

VIUTE(*,180)RCMATL,TOTTIN/24.0,TF)K},TPRHT 
WUTEC*.*)’  ’ 

HRITB(*.*)'  IRITIAL-TXMFERATURE-PmiLE' 

WUTEC*,*)’  ’ 

MtlTE(*,U9)(XXX(J,S),m(J,S)-273.1S, 

A  J-1, RIFTS) 

WITE(*,*)’  ’ 

NRITEC*,*)’  TOP-SURFACE-OONSTARTS’ 

HEITE(*,*)’  ’ 

HRITEC*,*)’  EMIS8  ALBEDO  ’ 

MUTE(  • ,  200  )EP3H ,  SMALU 
WITEI*,*)’  ’ 

WITE(*,*)’  IRPUT-UTER-SPECIFICAIIORS’ 

MaTE(*,*)’UTER  TBICERESS  VERT.  GRID  THERMAL  OIFF  BEAT  CORD’ 
MRITEC*,*)’  RO.  CM  SPACE-CM  CM**2/HIH’ 

A  CAL/MIR-CM-X’ 

DO  099SO  J4-1,RCMATL 
MRITE(*,210)J4,TBK(J*).SFRQ(J*), 

A  ALPB(J4),FE(J4) 

00036  CCRTIRUE 

HRITEC*, •)’  ’ 

MtITE(*,*)’  IRPDT  BOTTOM  BOURDART  DATA’ 

WITE(*,*)'  ' 

IP(.ROI.(LFLUXr.EQ.O))  GO  TO  08038 
NRXTE(*,80)LFLn(r 
HUTB(* .  02)DPHHD-273 . 13 
t«lTE(*,*)’  • 

»«ITE(*,*)'  ’ 
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00  TO  8MSS 

99836  irC.Kn.dAtncy.LT.O))  00  TO  99937 
HUTI(*.90)LrUOCy 
MUTS(* ,  93)DnMl*697 . 9 
WUtB(*.«)'  * 

00  TO  99938 

99937  MaTB(*,90)Ln.iar 

MUTBC  • ,  93  )Dn9<l*e97 . 6 

MUn(*,*)' - BOTTOM  SOtrACB -  SOBFilCE  BBHZATB  ATSSPMX  TIMP' 

MUTB(*,*)'EMISSmrT  GBOM  SBAFE  EHISSIVITY  GEOM  SHAPE  DEO  C' 
HRITEC*,*)'  PACT<0.-1.)  PACT(0.-1.)' 

MUTB(*,97)aEP,BE,ltZP,RE,TR-E73.13 
WaTE<»,»)'  • 

99939  tIUTEC*,*)'  ‘ 

IP(IVBG.EQ.O)  00  TO  190986 
MUTE(*,*)'  TEOETATIOH  PARAICTEBS' 

HEITE<*.»)‘  ■ 

mine,*)'  ooveraoe  state  oassmir  albedo 

A  POLIAOE  BEIGn* 

WITE<*.*)'  <0.0  -1.0)  <0.0  -1.0)  <0. 0-1.0) 

A  <CM)' 

l•UTE<* .  240  )SIiar ,  STATE ,  EPr ,  PQLA,  DOL 
W1TB<*.*)‘  ' 

)«TB<*.*)*  • 

»omE<*.*)’  • 

00  TO  198986 

C - 

99983  (XBiriHUE 

C  TO  CALCULATE-IHSOLATICM-OB-SLOPE-SOBFACE 

C 

C 

C  SOLPE-SOLAR-EEBITH 

ASSIGH  99931  TO  109932 
00  TO  08932 

C 

C  SOLVE-SOLAK-AZIMDTB 
00031  ASSIOK  09040  TO  100030 
00  TO  90030 

c 

C  CALCULATE-SLOPB-ATN08-ATTEll-AHD-<3.OUD-ADJUSTHEIITS 
09949  ASSIOH  00947  TO  X99946 
00  TO  09948 
C 

09947  OORTIBUE 
C 

98933  00  TO  109983 

C - 

09833  OOMTIHUE 
C  TO  ZEBO-TAUABUS 

c 

X-0 

00  TO  X99B33 


C . - - - - - 

9SSS2  OOKtim 

C  TO  SOLVE-StHitt-ZniTB 

C 

TT(»^MC»(TO«,24.0) 

TO-2 . o*Pi* (DAT- 1 . 0 ) /aes . 0 

DECL-0 . 00681S-0 . 399912*008 (TO l-K) . 0702S7*SIH (TO ) 

A  -0 . 008739*006(2 . 0*T0  )*0 . 000907*SIH(2 . 0*T0  ) 

A  -0 . 002697*008 (3.0*TO)>0.001*80*Sia(3.0*TO) 

ELr-(LAT/iaO*PI) 

TOCR- (  TTIC/ 12*PI  I'f  PI 

IF  ( TIM31 .  GT .  2 .  *PI  )Tliait-TDCIl-2 .  *PI 

AA^X)8(DECL)*C06(ELP)*C06(TIMER) 

a»*8IH(DICL)*3IH(ELF) 

C-AA+BB 

2'^ACaS(C) 

00  TO  199932 

C - - - 

99930  OCMTimiE 

C  TO  a(X.VE-aaLAil-AZIMIITH 

c 

xiRM*>co8(oaa)*stH(Ti»a) 

XDIO*<08(EU')*8III(DKL)-SIII(ELP)*C0S(TIMn) 

SAZ-ATAa(XinM/XDIKM) 

IF(.IKn.(XinM.LT.O.O.AIID.XDMCM.GT.O.O)>  00  TO  99944 

BAS-SAZ+PI 

00  TO  99943 

99944  IF(.IOT.(X1IUM.OT.O.O.AIID.XDIIOH.OT.O.O))  00  TO  99943 
SAZ-SAZ-PI 
99943  OORTimiE 
99943  00  TO  199930 

C - 

99948  (XUrrillUE 

C  TO  CALCOLAIE-SLOPE-ATMOS-ATTIM-AIID-CLOOD-ADJUSTMEIITS 
C 

SICF-C08(Z)*Oa8(8U)FE)*8IK(Z)*SIH(SL0FE) 

A  *C08(SAZ-SUItFAC) 

IF(.inT.(SICF.LT.O.O.OR.COS(Z).LE.O.O))  GO  TO  99941 
BOIHO.O 
00  TO  99942 
99941  l^l/008(Z) 

IF(.MOr.(M.n.O.O})  00  ID  99939 
IAL*0. 02023 

IF(DAT. 08.92.0  .AMD.  DAT. U.  132.0 )TAL— 0.02280 
11^3332 . 2/ (21 . 4-AU)6(ItB*ESAT(TA)  )  ) 

HATEi(-EXP(0 . 07074*(TD-273 . 13)4TAL) 

AB-0 . 271*(HATn*M)**0 . 303 

AO-0 . 083-0 .24 7*AL0010(FUS8/1000 .  *1 .  /M) 

AII81*(  ( 1 . -AB)*0 . 349*(  1 . -AO  )  /  ( 1 . -A0*0 . 2)*0 . 631 ) 

00  TO  99940 

99939  Al»l>1.0 

99940  Q^2.0*ANai 


QO  TO  99928 
99927  CORTIHUE 
99929  00  TO  199930 

C . - - - 

99981  COHTimiE 

C  TO  CALCULATE-SURFACE-AHD-LAm-TSIFEItATURE 
C 

C  SET-UP-IHITIAL-COHDITIOIIS 
ASSIGH  99921  TO  199922 
00  TO  99922 
C 

C  FRIRT-OUTFUT-HEAOIHG 

99921  ASSIGN  99919  TO  199920 
00  TO  99920 

C 

99919  OOHTIMUE 
C 

C  RUN-BEAT-FLOM-nOGSAM 
ASSIGH  99917  TO  199918 
GO  TO  99918 
C 

C  SET-UP-AHD-PRIHT-OUrm 
99917  ASSIGH  9991S  TO  199916 
GO  TO  99918 
C 

99915  IF(IRETRH.EQ.l)  GO  TO  4100 
GO  TO  99919 
4100  CORTIHUE 

GO  TO  199981 

C . . 

99922  CORTIHUE 

C  TO  SET-UP-IHITIAL-COHDITIOHS 
C 

PTTME-TOTTIM-24.0 

TINE-XXX(1,1) 

OIST-0. 

IFLAG-0 

IF  (TFHQ.LE.O)  TFRC^-TOTTIM 
DELT-TFRQ/60 . 

ITIHE*tMXl  (T0TTIM/DELT4 .9,1.1) 

IX-1 

IY-1 

GO  TO  99913 

99914  IF(IX.GT.NOMATL)  GO  TO  99912 
99913  IHTR(IX)-iy 

IF  <SFRQ(IX).LE.O.)  SFRq(IX)-TBK(IX)/10. 
HX(IX)^tAXl(THX(IX)/SFRQ(IX)4.9.1.1) 

RR( IX )-eO . 0*DELT/ (SFRQ( IX } •SFRQ< IX) ) 

C 

C  CALCULATICR  OF  HUSXICAI.  STABILITY  PARAMETER,  TO 
C  ENSURE  HUMERICAL  STABILITY. 

C 


STAn(IX)>AI.FH(IX)*i«(IX) 
ir(STABil(IX)  .81. 0 .  S)THai 
STABII(IX)-O.S 
ERD  IF 
C 

UYERS-O 
00  TO  9M10 

99911  tFCLATERS.GT.WdX))  00  TO  99909 
99910  XR-OIST 

RTABL-3 

OEPTa(iy)-» 

c 

C  m-TABLE-VALUES 

ASSIGN  99908  TO  199930 
00  TO  99930 
C 

09908  TEMP-TN 

STCR(l.iy)-TEMP 

8T0it(S,iy)-TEMF 

ST0it(6,iy)-FK(IX) 

8T0R(7,iy)-4IB0C(IX) 

STOR<4,iy)-0. 

3T0it(2,iy)-PK(IX) 

ST0il(3,iy)-ItB0C(IX) 

iy-iy+1 

OIST-OIST-I^FIiQ(IX} 
LAyERS-LAyESS^l 
00  TO  99911 
90909  IX-IX-M 

DIST-OIST-SFIKX IX- 1 ) 

00  TO  09914 

99912  JMAX>IY-1 
IllTR(IX)•^mAX 

■FIWT^IAXl  (TFIWT/TFIUH .  9 , 1 . 1 ) 

iJwrr-RFRiiT 

00  TO  199922 

c- - - 

09920  CXXrtlHUE 

C  TO  FRIHT-OUTFUT-BEADIIIG 

C 

C  IF(0UTCD.EQ.1)00  TO  1610 

C  IF(IVEO.GT.O)  GO  TO  1420 

C  MtITE(*,3S0) 

C  NRITE(*,3e0} 

C  MRITE(*,370) 

C  1«1TE(*,380) 

C  00  TO  109920 
C1420  i«ITE(*,310) 

C  MtITE(*,320) 

C  HUTE(*.330) 

C  MtITE(*,340) 

C  00  TO  109920 


w 


C1610  MUTE(*,390) 

SO  TO  I8SS20 

C . — - - 

99918  COMTimiE 

C  TO  Rmt-HEAT-FLOH-FROGiUM 

C 

IF{.IOT.(inJUS2.EQ.O))  SO  TO  99907 
ITER-HIT 
TIME-TIME+DELT 
99907  ZZA-ST0R(S,1) 

2ZB-ST0R(S,AMX) 

TEML-ZZA 

TDft-ZZB 

C 

C 

C  CALCULATE-BOUHDARy-OORDITIORS 

C  CALCULATE  EHERGY-BUDGET  TERHS  WITH  VEGETAIIOH  OR  WITHOUT 
C 

IF(IVEG.EQ.0)O0  TO  930 
ASSIGM  99903  TO  199800 
GO  TO  99800 

930  ASSIGR  99903  TO  199906 
GO  TO  99906 
C 

C  CALCULATE'UPFER-BOUNDARY-VALUES 

C  USING  ENERGY-BUDGET,  FIND  TOP  TEMPERATURE  WITH  VEGETATION 
C  OR  WITHOUT. 

C 

99903  IF(IVEG.EQ.O)  GO  TO  900 
assign  99903  TO  199797 
GO  TO  99797 

900  ASSIGN  99903  TO  I9990« 

GO  TO  99904 
C 

99903  IX-1 
J-2 

IMATL-HQMATL 

IF(.NOT.(NGMATL.BQ.l))  GO  TO  99801 
IZ-RX(IX)-1 

IF(.NOT.aZ.GT.O>}  GO  TO  99900 
C 

C  CALCULATE- INSIDE-MATERIAL- VALUES 

ASSIGN  99898  TO  199899 
GO  TO  89899 

C 

99698  CONTINUE 
89900  GO  TO  99902 
99901  GO  TO  99896 
99897  IF(IMATL.LT.l)  GO  TO  99893 
89896  IF(.NOT.(IMATL.GT.l))  GO  TO  99893 
IZ-HX(IX)-1 

IFC.NOT.dZ.GT.O))  GO  TO  99892 


c 

C  CALCULAR-IRSIDE  IMTERIAL-VALUE8 

ASSIOH  9989X  TO  I99S99 
00  TO  99899 

C 

99891  OCMTIIIUE 
C 

C  CALCUUTE-mnFiyCB-TALOES 

99892  ASSIGN  99889  TO  199890 
00  TO  99890 

99889  00  TO  99894 

99893  I>NX(Q()-1 
IF(.IIQT.(IZ.GT.O))  00  TO  99886 

C 

c  calculatb-ihside-haterial-values 

ASSIOH  99887  TO  199899 

c 

00  TO  99899 
09887  CONTINUE 
99888  CONTINUE 

99894  IHATL-IMAT1.-1 
00  TO  09897 

9980S  CONTINUE 
C 

C  CALCUUTE-LGNER-BOUHIMRy-VALUES 

09902  ASSIGN  99883  TO  199886 
00  TO  99886 

C 

90883  00  TO  199918 

C - - - - 

90906  CONTINUE 

C  TO  CALCULATE-BOUNDART-CONDITIONS 
B  -  -«(1) 

T-TIME 

C 

C  CALCULATE-BOTTOH-BOUHDART-BEAT-TEEte-AFiai-DFRM-BFRM 
ASSIGN  99880  TO  199881 
00  TO  09881 
C 
C 

C  ATMOSFBERIC-IHFIIASED-EHISSION-ATZIIH 
09880  ASSIGN  90878  TO  109879 
00  TO  90879 
C 

00878  CONTINUE 

C  CALCULATE-SOLAR-IHSaUTION-FOR-DAy-AND-TItE 
DAT-XXX(1,8) 

ASSIGN  42  TO  100069 
00  TO  09989 
42  BTB1B8  8UN 

IT  (BTERM.OT .  0)BTEi9t*BTEBM*9HALLA 


C 
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C  CALCULATE-OaMVECTIOH-BTOIM 
itSSIOR  88876  10  188877 
00  TO  88877 
C 
C 

C  CMLCULATE-EVAFCRATIVE-HEAI-LaSS-DTEKM 

88878  ASSIGN  88874  TO  188873 
00  TO  88873 
C 

89874  D  -  ATERH  -4  BTESM  -  BTEm-DTEBM 
321  CaMTIHUE 

00  TO  188808 

C - - - - - 

99881  oammiE 

C  TO  CALCULATE-BOTTOH-BOOIlDAXy-HEAT-TESMS-AniH-DnM-Bm 
C 

BFI»»>TB 

IF(.ini.(TB.BQ.O.O))  00  TO  99872 

AFIB^l.O 

DFBH-DFBMO 

00  TO  98873 

89872  AniH*FACTE*TEMt*TD«*TEMl 
DFMt-OFmi 
99873  00  TO  199881 


99879  OOMTIHIIB 

C  TO  ATMOSFHERIC-IMFRARED-ZMISSIOH-ATQiM 
XH-TIME 
IITABL-2 

C 

C  (an-TABLE-TALDCS 

ASSIGN  88887  TO  188830 
GO  TO  98830 

C 

89867  XB-TH 

XH-TIHE 

NTABL-10 

C 

C  GET-TABU-VALUES 

ASSIGN  99866  TO  199930 
00  TO  99930 

C 

99866  TA-TN 

XH>TIME 

IITABL-3 

C 

C  GET-TABU'TALUES 

ASSIGN  99863  TO  199930 
00  TO  88930 

C 

88863  CLOUD-TN 
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I«>(TAK-273.13) 

SA-e.  108*lia*EXP(  (AC*TAC)/  (TAK-BC) ) 

ALFBI»(0.61>0 .  OS*SQRT(EA)  )*(  1 . 0-f  (CU(HCU»D)*(CL0UD**2) ) ) 
DONMn-0 . 8132Z-  10*TAX**4*ALFHI 
ATERH-OOHHIR 
00  TO  199879 

C - - - - 

99877  OCRTIIIOB 

C  TO  CALCUUTE-COMVECTIOM-IITESM 
C 

XM-TIME 

RTABL-e 

C 

C  GET-TABU-VALUES 

ASSIGH  99862  TO  199930 
GO  TO  99930 
C 

99862  SFEEO-YII 
TAK-TA 
ZASB-ZA 
TSX-TEHL 

RBDA— 0 . 001*0 . 348'FRESS/TAIC 
1200  THETAZ-TAK*FACTB 
THETAS-tSK*FACTH 
DTBETA-(TBETAZ-TBETAS)/ZASa 
DU-SFEED/ZASB 
TBETAV- (TBETAZ*TBETAS ) /2 . 0 
RI-0*DIHETA/ (THETAV«OU**2) 

COEWS.O 

COE2-1.175 

EX-.75 

IF(TSE.GT.TAX)00  TO  31 

IF  (RI .  GT .  0 . 2  )IU> .  19999 

COTl-5.0 

COE2-1.0 

EX'>2.0 

31  BTER-ilB0A*KSQ*ZASB**2*DU 
&  *{COE2»(1.0-OOE1*IU)«*EX) 

BTEIBMTER*CP*DTHETA 
99864  GO  TO  199877 

C . - - - 

9987S  COMTIHUE 

C  TO  CALCULATE-EVAPCRATIVE-HEAT-LaSS-DTEIH 
C 

XH^IHE 

HTABL-Z 

C 

c  OT-TABU-VALUES 

ASSIGN  99839  TO  199930 
00  TO  99930 
C 

99839  RI-TN 
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XM-TOC 

KTABL-IO 

C 

C  aET-TAXLE-VALOES 

ASSIGM  9MM  TO  199930 
00  TO  99930 

C 

99838  CTEMA-n 

KIDIPA-CTDM 

CTEHA-CTEMA-273.1S 

KTDffOTOO. 

ES-EEP(  (AC*(KTEHFG-273 . 13 ) )  /  (KTSCG-BC)  )«6 . 1071 

EA>EEP(  (AC*CTEMA)  /  <KTEMPA-BC)  )*6 . 1071>SB 

DO-0 . 622/FRESS* (EA-ES ) *UET/ ZA 

XL-397 . 3-0 . 386*(CTEMA-HCTEHFG-273 . 13)  /2 . 0 

DTEIIM*BTER*XL*D6 

00  TO  99061 

99860  DTXHt^O.O 

99861  GO  TO  I9987S 

C . . . - . . . - . 

99904  CORIIIUE 

C  to  CAUUUTE-OmR-BOUimARY-VALUES 

Tl-8TABI»(l)*(STC»<l,3)-2.»STC*<l,2)+STCIl(l,l))+STQR(l,2> 

III-O 

830  lIl-m+1 

T2-STaR(  3 , 1  )****FACTA*SPIK}(  1  )+FK<  1  )*ST0R(  3,1) 

A  -<FK(1)»T1+D*SFRI}<1)) 

T2-T2/ <  * .  •FACTA»SFRI)<  1  )*ST«l(  3 . 1  )»*3+nC(  1 ) -SFR0(  1 )  *DDDT ) 

STCR(3,1)-ST0R(3,1)-T2 

TEML<«TaR(3,l) 

C 

C  SURFACE  iTORE  (TBML,  ST(Xt(3.1))  SOI  ALLOWED  TO  RISE 

C  ABOVE  0  3T  SHOW  PRESEHT.  THE  SHOWPACK  IS  SATURATED 

C  AT  O'C  AMD  UHSATURATED  AT  TEMPERATURES  BELOW  O'C. 

C 

IF(T»L.GE.273.13)THER 

TIML-273.13 

WBT-1.0 

ELSE 

WET-0.0 

END  IF 

IF(STQR(  3 , 1 )  .6T .  273 . 13)TBE1I 
ST0R(3,1)-273.13 
EBD  IF 
C 

6TEIB*>-FK(  1  )*  (ST0R(  3 , 1 ) -T1 ) /SFRQ(  1 ) 

ASSIGR  823  ID  199877 
GO  TO  89077 

823  ASSIGII  810  ID  199873 
00  TO  99873 

810  DOEWi^TEWHBTBRM-HlEIIM-UlEIIM 

tF(AB8(T2).LT. 0.003  .OR.  III. GT. 30)00  TO  199904 
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i 

I 


ODOT— (nm-om 
MMBI 
00  10  630 
C 

C - 

99880  OQRTimiE 

C  TO  CALCUlAIE-llBIOE-HATEmL-VALUES 
C 

00  TO  908S6 

998S7  IFdZ.LB.O)  00  TO  098SS 

00836  conim 

STaR(3.J)>4Tat(l.J)>STABII(IX)*(STCR(l.J'l)-2.*STn(l,J) 

8  >aToo(i,jn)) 

J-J+1 
IZ-IZ-1 
00  TO  09837 
99833  00  TO  190890 

C - - - - 

C  TO  OtLCULATE-IirrEllFACE-VALUES 
C 

99890  OORTIHin! 

■OOC9^TOR<  6 ,  J- 1 )  /  8FBQ  ( IX  ) 

DCtXF^TORC  6 ,  J-M  ) /SF1IQ(  IXn ) 

OOaeP-BCQEF-HXXKF 

ACOeP-BCOEF/ (2.  *STABII(1X)  )'H)aSF/ (2 .  •STABHCmi)  > 
STOR(3,J)-ST(*{l,J)+(BCOTF*STOR<l,J-l)-CCOEr*STOR{l,J)+DCOEF« 
&  STCR(1,J4^2})/«CXSP 
ST0R(  3 ,  J-M  )-3T0R(  3 ,  J ) 

IX'IX-M 

J-J+2 

00  TO  199800 

C - - - 

C  TO  atLCULAIX-LOHER-BOOinMRr'VALaES 
99886  ir(Ln.UXT.EQ.O)  00  TO  880 
I-l 

X2-FACTD 
870  conmuE 

Rl-SiaH8*BEF*K*8T0R(3,J)**« 

01— FK(IICmTL)*(8TaR(3,J)-STCR(l.J-l))/SFRQ(lRH8TL) 
F2-*.0*SIOM8*BEP«BX«ST(Il<3,J)»*3-FK{BOM8TL)/SFRQ(HOM8TL) 

F2f  -  (BZ-Rl-tOl-fDFMO/FZ 
ST0R(3,J)-ST0R(3,J)  *  F2 
I-I+l 

IF(ABS(F2).6T.0.01  .MD.  I.LE.30)  00  TO  870 
880  IF(LFLtncy.n).O)  STai(3,J)^TCR(3,J) 

00  TO  199886 

C - 

90916  CORIiniE 

C  TO  SBT-DF-AUD-fMinPUUTE  OUTPUT 
C 


IFL8O2-0 
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irCITER.LE.l)  GO  TO  124S 
ITBS-im-l 

IF  (lEFSHT.R.O)  00  TO  1328 
IF  (IFnr.LE.l)  00  TO  12*4 
IF  (ITDC.GT.l)  00  TO  1328 
12*4  CCHTimiE 
C 

C  HARIFULATE  OOTFUT 

ASSISI  89846  TO  199847 
00  TO  99847 
C 

99846  00  TO  1328 

1243  IF  (ITnC.GT.l)  00  TO  1269 
C 

C  MAlllFULATE  OOFUT 

ASSIOH  99843  TO  199847 
00  TO  69847 
C 

99843  niETIW-l 
00  TO  1333 

1289  ITDB-lTIME-l 

IF  (IFRin.LB.l)  00  TO  1279 
inUrflFiWT-l 
00  TO  1303 
1279  CORTIiniE 
C 

C  MAHIFDLATE  OOTFUT 

ASSIOH  99844  TO  199847 
00  TO  99647 
C 

99844  IFRHT-HFiUlT 

1303  J>1 

IZ-HOMATL 
1306  nC-HXdD+l 
1311  STCR(l,J)>STai(3,J) 

ST0ll(2.J)-ST0il(6,J) 

ST0R(3,J)>ST0R(7,J) 

J-J+1 

IF  (IX.U.l)  00  TO  1329 

Dt-nc-i 

00  TO  1311 

1329  IF(IZ.LB.l)  00  TO  1333 
IZ  -  lZ-1 
00  TO  1306 
1328  IFU62-1 
1333  COHTIMUE 

00  TO  199918 

C - - 

99847  OOHTimiE 

C  OUTFOT  MAHIFOLATICH 

C  BOURLT  BnSGr  VALDES  FOR  BOTH  UHVERTATZD  AHD  OHVERTATED  SUHFACES 
C  ARE  OORVERTED  INTO  lyllE*^.  THE  EREROT  VALUES,  SURFACE  AMD  BASAL 
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C  TEMPBBATUUS  MD  COHUSPONDIlia  KXIR  OF  TBE  DAY  ASZ  PLACED  IH  A  DATA 
C  HAim.  DATW. 

C 

TPACK-1.0 

IF(.im.(TIIC.GT.PTYIB))  GO  TO  99843 
DO  99842  JiaC>l,KHATL41 

TITLE(JEX)-(ST(m(S.IJ)-273. 13) 

99842  ooimm 

IIDX>TIia 

IF(ian(.BQ.O)IIDX-l 

c 

C  RESULTS  CALCOUTED  M  ealea'^Bln'^  IE.  lyiBlo'^  AMD 
C  CmiVBKTBD  TO  lyhx'^  BY  *80.  THE  DAILY  TOTAL  IS  THEM  DIVIDED 
C  BY  24  (ELSEMBERE  IR  THE  PROGRAM)  TO  HAVE  AHSMER  IM 
C  IS  MM  THE  HDRE  WIDELY  USED  UHIT  ARD  THEREFORE  IS 

C  USED. 

C 

STP-TITU(l) 

BOTTP-STOR( 3 , J ) -273 . 15 
K10HC1041 

IP(IVES.Bq.l)  00  TO  1110 
C 

C  ERERGY  VALUES  FOR  URVEGETATED  SURFACE,  COHVERTED  TO  lyhr'^. 

C 

IGBR^(0 . 813E-10*EPSR*STCR( 3 , 1 )**4 )*60 

I9QL-<BTERH/aMALLA)*60 

IABSCR-(ISOL*SMALLA) 

IATERI^ATERM*eO 

IBTERN4TERM*60 

IDTEi»^DTERM*80 

C 

1110  CORTIRUE 

DAIMECKIO, 1)-TIHE 
DAIHX(K10,2)-TirLE(l) 

C 

C  CALCULATE  VE6ETATICR  RADIARCE  VALUES  IF  VEGETATIOR  MODELLED 
C 

ASSIGR  1400  TO  11410 
IF(IVBG.EQ.l)  GO  TO  1410 
1400  CORTIRUE 
C 

IFIIVEG.EQ.  OTHER 

C*********************************************************************** 

C  OPTIOR  2-BOORLY  VALUES  FOR  ERERGY-BUDGET  VALUES. 

C  IRSERT: 

MRITE(9, 1BO)AMOO(TIME,24 . ) ,TITLE(1) , IBGR, ISOL, lABSOR, 

A  lATERM.IBTERM.IDTERN 

C 

DATMC(K10,3)-iaBR 

DATMC(E10,4)>I80L 
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MIHC(U0,3>-UaSGB 
DAm((Kio,e)-iATnM 
DAK«(K10,7)*IBniM 
DAIMKKIO .  a>>IDniM 
DA1M((K10,9)>BOTTP 
ELSE 

C***********«******#*«****************************************i 

C  OPTIOH  2-aoa8LY  VALUES  FOR  BHEBSy-BUDGET  VALUES. 

C  INSERT; 

HR1TE(9 , 270 )AMX(TIME . 24 .  ) , ISURFtHXREFRA,  ISURFG, 

A  TEFFR-273 . IS , TEFP-273 . 13 , TlML-273 .15,Tr-273.15,ISOL.SG 

C*************************************************************. 

c 

0A1ME(E10 , 3  )««niDaB*«0 
DAlt«(E10,4)-SG 
0ATMX(K10, S)-IAS80K 
DATI«(K10 , 6  )<«LD*eO 
DAna(K10,7)-HS6*e0 
OATM((K10, 8)-ELG*XLl*60 
DAn«(K10,9)«BOnP 
END  IP 

99843  00  TO  199847 

C . . . . 

99800  COHTINUE 

C  TO  CALCULATE-BOUNDARy-COiailTIOIIS-WITB-VEG 
T-TIW 
XI»>riME 
IITABL>e 
c 

C  ^-TABLE-VALUES 

ASSIGN  960  TO  199930 
00  TO  99930 
C 

960  UA>TN 

xH^na 

IITABL-4 

c 

C  AIMOSPBERIC-INPRARED-IMISSIOH-ATERH 

ASSIGN  980  TO  199879 
00  TO  99879 
C 

C  CALCULATE  INCOMING  SOLAR  RADIATION. 

960  CONTINUE 

0AT4CXX(1,8) 

ASSIGN  43  TO  199983 
GO  TO  99983 
43  SOL-SUN 

C 

IP(UA.LT.10.0)UA-10.0 

UAF-0. 83*SiaPMIA*S9tT(CIIB)4(  1 . -SI(9)*UA 

OELIMP-3. 

CP-0.01»(1.+30.0/UAP) 
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Du-(iM-aAr)/u 

■8-1/  ( .  09+ .  oo2i*<aat*u ) ) 

■C-IS+STATB/ ( 7 . o-sior > 

Air(i)-Tr 

ASSIOa  1210  TO  1090 
00  TO  090 
1210  COHTim 

nB(i)-rm 

■DDC-O 

1240  TF-TF+DELlMP 

■DEX-aOEX+l 
ASSia  1220  TO  1090 
00  TO  090 
1220  CCtniMUE 

m  RDEX .  EQ .  3 .  AMD .  mo .  LT .  0 . 0  )  THEM 

mo-mo*(-i.o) 

Em  IP 
PEB(2)-FEIO 

IP(FEB(l)+FEa(2).LT.0.0)  00  TO  1230 

IP(1IDEE.EQ.3)TBER 

OOTO  1230 

EHD  IP 

IP(ABS(PEB(2) )  .or .  ABS(PEB(  1 ) )  IDELTMP— 9 . 

IPdOEX.LT.  100)00  TO  1240 

l«irE(*,*)’P0LIA<O  EHEROT  BUDGET  HAS  HOT  C3IOSSED  X-AXIS* 
NRITEC+.-l'APTER  100  SEARCH  STEFS.  CHECK  IIIFUT  DATA.' 
STOP 

1230  CQHTIIIUE 

RDEXl-O 
ATF<2)-TF 

1270  SLOFE2-(FEB(2)-FEB(l))/(ATP(2)-ATF(l)) 

BIirr-PEB(l)-SLOPE2*ATF(  1) 

TPO— BIIIT/SL0FE2 
HDEXl-HDEXl+l 

IF(ABS(TP-TFO).LE. 0.01)00  TO  1260 

IP(IIDEX1.EQ.3)TBEH 

OOTO  1260 

EHD  IP 

TP-TPO 

ASSIGH  1290  TO  1890 
00  TO  990 
1290  OOMTIHUE 

IF(PEIIB*PEB(2)  .or .  0 . 0 )  IP-2 
IP(PEIIB*PEB(1)  .GT.0. 0)1^1 
ATP(IP>-IF 
PXB(IP)-PE1IB 
00  TO  1270 

1260  00  TO  I99B00 

C - - - - - 

C  TO  CALCUUTE-aPPSR-BOtmDAinr-VALOES-FCR-PXIAGB 
89797  OOllTIIIUE 

ISLIIIF-9. 
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1310 


1340 


1320 


1330 


1370 


1350 


1360 


C— - 
C 

050 


Air(l)>TBC 

ASSICM  1310  TO  11300 

00  TO  1300 

CORIHOI 

rSBC  1)400 

100(4 

TDf.-Ta«.-fOQ.1MP 

IOD(HOO(41 

ASSIGN  1320  TO  11300 

00  TO  1300 

(XMTIinK 

FEB(2)-FOa 

IF(m(l)*m(2).LT.0.0)  00  TO  1330 
IP(RID(.n).S)TBDI 
OOTO  1330 
OO  IP 

IP(ABS(PEB(2)).eT.ABS(PEB(l)))DEL1MP— 5. 

IPdOEX.LT.  100)00  TO  1340 

HUTe(*,*)‘aBOUD  DIOIST  BUDGET  BAS  HOT  CROSSED  X-AXIS’ 

MRITE(*,*)‘APTER  100  SEARCH  STEPS.  CHECK  IHFUT  DATA.’ 

STOP 

CONTimiE 

HDEX14 

ATP(2)-TDC. 

SU)PE2-(PEB(2)-PEB(1))/(AIF(2)-ATF(1)) 

BIHT>PEB  ( 1 } -SLOPBRMTP  ( 1 ) 

TPO— BIIIT/SU)PE2 
HDEXl-HDEXl-fl 

IP(ABS(TBM.-TPO).LE. 0.001)00  TO  1360 

IP(H0EX1.BQ.3)TBEH 

GOTO  1360 

EHD  IP 

TBa.«TP0 

ASSIGN  1350  TO  11300 
00  TO  1300 
CONTINUE 

IP(PEHB*PEB(2) .OT.O.O)IP-2 
IP(PENB*PEB( 1 ) .or . 0 . 0 )IP-1 
ATP(IP)-TDa. 

PEB(IP)>PEHB 

00  TO  1370 

STat(5,l)>TB1L 

lP(TOa..aE.273. 13)THEH 

IIHL-273.13 

EHD  IP 

00  TO  IS0797 


TO  CALCULATB-EHERGT-BUDGET 

TAP- ( 1 . -Slop )*TA4SI0F* ( 0 . 3-TA+O . 6*IP40 . 1-TEML ) 
0T8BTA-(TA-17)*PACTH/XA 
THETAV-(TA+TP)*PACTH/2 . 0 
RI4*DmTA/  (THETAV»DO«2  ) 
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aaCMT—0 . 001* .  346*mESS/  ( (Tr*IA)  /2 .  ) 

OQBl'13. 

OaK2>1.17S 

BC-.7S 

ir(si.u.o.)ao  to  uso 

IF(U.GT.0.2)U>KI.10« 

COBl-S 

0CB2-1. 

SX-2.0 

1280  ctmiKia 

BTEB*«aaAr*KSQ*ZA**2*DU 
&  •COe2*(l.-OOEl*IU)**EX 
BSF-aTER*CP*DTBRA*60 . 

XL-307.3-0.36«*rAr 

BA-(AL06(  (2A-Z08P) /20  )*C0E2*(  ( 1 . -C0E1*RI  )**EX)  )**2 
&  /(.18*UA) 

SOP-HA/ (BS-fRA) 

OFSIDPmSAT<TF)+<  1 . -IU)P)*<JAP 

f'AF-<  1 .  -SI6F)*q(TA)*SI(7*(Q(TA)*0 . 3*<3P*0 . 6-K}G*0 . 1  > 

JF— (RHOAF*CP/0.66)*(ESAT(TP)-E(TA))/{RA+1IC)*60. 

IPCEF.LI.0.0)EF-0.0 

SaM^FOLA^SOL. 

XU»»'EPF*ATE8M 

TG4-EPF*EPSa/EPl*SIGMA*TEML*** 

TF*-(EPH-EPS* )  /EP1*EPF*SICMA*TF**4 
FEIIB-SIW*<SaiBHiajlC»HTG*-TF*)-HSr-EF 
C»  TO  I9S0 

C - 

C  TO  CALCOLATE'EREilST-BlIDGET-FCR-GnOUIID 

1300  OOMTimn! 

T1-ALPB(  1  )•]«(  1  )*(ST(ll(  1 . 3 )'2 . *STOR(l . 2)*STCXt(  1 . 1)  ) 

&  +STCXl(l,2) 

IF*-SIGHA*TF*** 

TG4-SIGHA*TEHL*** 

QI>4*eT*q8AT(TEML)+(l. -WET)*QAF 

SBCMB-0 . 001*0 . 3*8*PltESS/TAF 

XLl-597 . 3-0 . 568* (TAF*TEML-2 . 0*273 . 15) /2 . 

S®-<l.-SIGF)*SOt 

BLU-C 1 .  -SIW)*<EPSII*TG4+(  1 .  -EPSH  )*ATERM) 

&  +SIGF*(EPS»*IG*+(1.-EPS«)*EPF*TF*)/EP1 

RU)-(  1 .  -SIW)*ATEIIM*SI(ar*(EPF*TF4+<l .  -EPF)*EPSH*TG*  )  /EPl 
BS0-«BQA0*CP*CB6*UAF*(TEM.-TAF )*60 . 

ELO*R8aMC*CaS*(MF*  (Q6-QAF )  *60 . 

FEIIB-SMALLA*86-RLU*RLO-BSO-EL6*XL1*(T1'TEML)/SFRQ(  1)*FK(  1) 
00  TO  11300 

C - - - - - 

C  TO  CALCULATE-BADIARCE-TALUES 
1410  oonriiraE 

BZFBAIHC ( 1 . -8I0F)*( l-EPai)+810F*( 1-EPF) )*DOM»IP 
FO(ja»*EPF*0 . 8132B-  10*TF**4 
GRin)eB-EPSII*0 . 8132E-  10*TIK.**4 
S0BFa»4ieP*F0tjra*(  1 . -SIOF)*(aUIDGB 


EEr-SlGP«EPP+<  1 .  -SI(7)*EPSM 
TEFF-(SURFGB/0 . ei32E-10 )•• . 23 
ISURro-SURFG8*60 

TEm-(  (  (SURFGB-i-REFIiAO) / (0 .  S132E- 10  >  >** .  23 ) 

ISEFItA-ItEFIlAD*eO 

ISOL-SOL*60 

SO>(  ( 1 . 0-SIGT  )*SOL)*60 
IABSaR-SMAlXA*SG 
00  TO  lUlO 

c - 

EHD 


c*****************************************************************. 

C  OFTIOH  3-MASS  OOHSBIVATIOH  ROOTIIIE 
C 

C  FOR  A  2-LAYER  PACE 

3001  OEUtIIS-(DEIjqM/(LHEATF*VIIEHSM(l,6)*1000)))*iaO 
IP  (DELSRS .  GT .  0 . 0  )TBEH 

IFCDELRlIS .  (a .  TBKMa  ( 1 . 1)  )THEH 
ABC-OELSRS-TBICMa  ( 1 , 1 ) 

DFOTBEMXEd.l) 

THKHUd.D-O.O 

OCLqM>(  ABC*LBEATF*  (  VKEWSRd .  6  >*1000  )  )  )  / 100 
OELROS>(DELQM/(LBEATF*(VOLOSNd.6)*1000)))*100 
IF(DELR0S.GE.TB1CMX2(2, 1))THEN 
GOTO  3002 
ELSE 

TH1CM)C2(2. 1)-THKKIC2(2.  l)-DELROS 
DEPTB»THXK(2(2.1) 

HOMATL-l 
IHPIS-2 
DO  800  J-1,6 
THKKXK  I)*VOLOSHd,J) 

800  OQRTimiE 

TBKmid,l)*OEPTB 
DEFKXld,2)*AVSTF 
OEPHXl(2,2)>AVBOTTF 
DEmXl(2,l)*CEFTH 
GOTO  3000 
ERD  IF 
ELSE 

TBXHE2 (1 , 1 )*TBKKX2 (1 , 1 ) -DELRHS 

0EmE2d,2)-AVSTF 

DEFMX2(  3 , 2  )-AVBOTTF 

DEFH(2(3,1)~0EFTB 

OEiMC2  (2 , 1  )~TB1EMX2  d ,  1 ) 

021(02(2 , 2)-<  (DER02(2, 1)/DEPTH)*(AVB0TTP-AVSTP)  )+AVSTP 
GOTO  3000 
ERD  IF 
ELSE 

GOTO  3000 
ERD  IF 
C 

C  FOR  A  1-LAYER  PACK 

4000  IF(TIIKI01d,6).EQ.VREHSRd.6))TBER 

DELRRS-(DELQH/ <LBEATF*(VREHSRd.e)*1000>  )  >*100 
1F(DELSRS .OT . 0 . 0 ITBER 
IF(DELSRS .  OE .  THKIOl  (1 , 1 )  ITBER 
GOTO  3002 
ELSE 

TBBOl  (1 ,  D^TBKMCl  (1 . 1 ) -OELRRS 
DEFT»4B901d,l) 

ISf(01d,2>-AFSTF 
DBIIOl  (  2 , 2  )*ATBOITF 
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DER«1(2,  D-TBUHCKl,  1) 

GOTO  5000 
EMD  IP 
ELSE 

GOTO  5000 
EHD  IF 

DELBOS-(DEL(]H/(LBEATF*(VXS81I(1.6)*1000)>)*100 
IF(DELIIOS  .GT .  0 . 0  )TBER 

IF  (DELROS .  GE .  TBIEMXl  ( 1 , 1 )  )THEII 
GOTO  5002 

THKMXld,  D-TBEMXKl.D-DELSOS 
DEPTB-TBIKCKl.l) 

DEFHX1(1,2)>AV8TP 
DEFHX1(2,2)-AVB0TTP 
DEPHCl  ( 2 , 1  )-THXH(1  ( 1 , 1  > 

GOTO  5000 
EHD  IF 
ELSE 

GOTO  5000 
EHD  IF 
EHD  IF 
C 

5000  IF(HaMATL.BQ.2.AH0.TBXME2(l,l)-Lr.l.0}TBEH 

DEPTB-THKMC2  ( 1 , 1 ) +TBENU  ( 2 , 1 ) 

HOMATL-l 
HIPTS-2 
DO  A50  J-1,6 
THKHEl  (1 ,  J)-TOLDSH(  1 .  J) 

*50  aUTIHUE 

THEHXKl.D'DEPTB 
DEP»«1(2,1)"DEFTH 
DEmXl(l,2)-AVSTF 
DEP(«l(2.2)-AVBOTTF 
EHD  IF 
C 

MiSH-DPG+OELROS 
HEHS>DFO*VHEHSH  ( 1 , 6  } 

HE08M>ELitOS*VOLOSH  ( 1 ,  e  ) 

MiHE>WEHS*ilE06 
IF(t«SH.LT.0.0)THEH 
MISH^.O 
EHD  IF 

IF(»IIHE.LT.O.O)THEH 

MINE>0.0 

EHD  IF 

C**«***************************************************< 
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VARIABLE  OEFINITIOHS:  SKMO 

(1)  MAIN  PROGRAM 

AIRTP  Alt  tMparatur*,  *C. 

-2  -1 

AVABSR  Maan  dally  ahortwava  radiation  abaorbad  at  tha  aurlaca,  HIb  day 

-2  -1 

AVATER  Maan  dally  ahortwava  radiation  abaorbad  at  tha  auxfaca,  MJm  day 

AVBTTP  Naan  dally  taa^ratura  at  tha  baaa  of  tha  tnowpaek,  *C. 

_  -2  “1 

AVDTER  Maan  dally  avaporatlva  haat  flux.  Him  day 

_  -2  -1 

AVOR  Maan  daily  raflaetad  loogwava  radiation,  MJm  day 

AVHTER  Maan  daily  aanalbla  haat  flux,  HJa  ^day 

-2  *1 

AVSOL  Maan  dally  Ineoadnc  ahortwava  radiation,  tUa  day 

AVSTF  Maan  anow  ta^fMratura ,  *C. 

BOnP  Ta^Mratuxa  at  tha  baaa  of  tha  snowpack,  *C. 

CELL  Idantlficatlon  n«ad>or  of  tha  coa^tatlonal  coll  bains  Bodallad. 

CRITDF  Critical  anowdaptb.  Sea. 

DAIMX  Matrix  holding  TS1M  output  valuaa  for  aunipulatlon . 

DAY  Julian  data  of  day  balng  aodallad. 

DAYl  Julian  data  of  SIKMO  Initiation. 

DDATE  Julian  data  of  day  bains  aodallad. 

-2  -1 

DELQS  AQ  ,  anowpack  Intamal  anarsy  chansa,  MJb  day 

_  -3 

DEIISM  Danalty  of  anow,  $cm 

-3 

DEHH  Danalty  of  watar,  gem 

DEFHXl  Oaptb  natrlx  for  a  I'layar  pack. 

0EFMX2  Dopth  aatrlx  tor  a  2-layar  pack. 

DEPIB  Snowpaek  dapth,  oa. 

DEFTBl  Oaad  In  calculation  of  dapth  whan  convartlns  'naw  anow’  to  'old  anow’. 

dU/dt,  part  of  tha  anowpack  anarsy*budsat  aquation  (aqulvalant  to  dQ^),  radundant. 


DODT 
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EHAVH 

GTERH 

BCAFSN 

BCAPH 

ILAFSE 

IKAIH 

ISHOH 

IVEG 

K1 

K2 

K3 

X25 

K30 

KELVIH 

XaHTRL 

LBEATF 


-2  -1 

Enazsr  avalliabl*  for  salt,  calculatad  ualnc  anargy  totals,  MJa  day  *. 

*2  '1 

Ground  haat  flux,  calea  Bin 

Baat  capacity  of  snow,  kJkg  ^*C 
Baat  capacity  of  watac,  kJkg  ^*C 

OatazBinss  which  aiavation  data  fils  is  ussd  to  laodify  ths  sir  taopazstuces  with 

lapaa  zata. 

Indicatas  occuzanca  of  zain-on-snow.  If  izain^O  than  no  rain-on-snow  occuzad,  If 
izain-1  zain-on-anow  occuzad. 

Indicatas  if  snowad  on  tha  pravious  day,  isnow^O,  no  snowfall  on  previous  day, 
lsnow>0,  snowfall  on  pravious  day. 

DatszBinas  whsthar  VE6IE  is  activatsd  or  not  and  which  vagstation  input  fils  is 
ussd. 

Count  of  tha  nuBbaz  of  days  axidsllsd,  is.  ths  nund>sr  of  timss  tha  'do  5001'  loop  is 
ussd. 

X2~0  indicatas  that  SMGMO  aodsllad  soil  in  ths  pravious  Itszatlon  or  that  SHCMO  is 
not  aodalling  ths  first  iteration.  K2~l  indicatas  that  that  5KHO  modelled  snow  in 
ths  pravious  Itszation  or  that  it  is  tha  first  Iteration. 

K3~0,  indicatas  soil  sndallad  in  pravious  iteration,  IC3>1,  indicatas  snow  oodallsd 
in  previous  itszation. 

Count  of  the  days  slapssd  sines  ths  Isst  snowfall. 

Indicator  to  facilitate  ths  calculation  of  ssntp,  in  tha  swlt  calculations. 

Air  taaqpsratuza  eonvsrted  to  dsgraas  kslvln. 

Count  of  tha  nuid>ar  of  days  aa>dallad. 

Latent  heat  of  fusion,  HIkg 

Maltrata,  ea  snow.day  ^ 

Haltzata,  as  snow  wataz  squivalaot.day 


BBTIUD 


■at  radiation,  a  ^day 
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Him 

nCDIIT 

tavoL 

ppm 

PTERH 

SB 

SMDPl 

SNPPTH 

SHVOL 

SOILDP 

SOILTP 

SPBTI 

SP8TN 

SSOLTP 

STP 

TABSQR 

lATEBM 

TDTESH 

T6BR 

THKMXl 

THU«2 

THTESM 

TDCS 

KMX 

mix 

TSOL 


Th«  ntaibax  o(  tiaai  tha  output  tlao  print  froquoncr  !■  divldblo  by  th«  tliM  it«p«, 
usod  to  dotoxalna  «hon  to  print  output. 

Uaod  with  SDXTSOL  to  toxainoto  SXONO. 

Iho  niadbor  of  aatorial  layara  uaad  in  aolvins  tha  haat  flow. 

Praclpitation,  m  watar. 

-2  -1 

Baat  input  to  pack  by  rain-on-snow,  HIb  day 

Ralativa  hiaidity,  Z. 

Initial  anowdapth,  ea 
Praclpitation  aanunt,  cb  anow. 
araa  x  dapth,  B 

dapth  of  toil  profila  at  tha  baaa  of  tha  anowpack,  cb. 

Soil  taaparatuxa  at  dapth,  takan  at  aoildp,  ’C. 

Spacific  haat  of  ica,  2.10x10^  Jkg 
Spaciflc  haat  of  watar,  t. 18x10^  Jkg  ^*K 
Suxfaca  aoll  taq^ratura,  *C. 

Surfaca  taaggMxatuxa,  *C. 

-2  -1 

Total  daily  shortwava  radiation  absorbad  at  tha  aurfaca,  HJb  day 

-2  -1 

Total  daily  incOBlng  longwave  radiation,  MJb  day 

-2  -1 

Total  daily  avaporatlva  haat  flux,  MJb  day  . 

Total  dally  raflactad  radiation,  HIb  ^day 
Thlcknaaa  aatrlx  for  a  1-layar  pack. 

Thicknaaa  aatrix  for  a  2-layar  pack. 

Total  daily  aanaibla  haat  flux,  HJm  ^day 

Dally  obaaxvatlon  tlaa,  uaad  to  aid  tha  aaquantlal  reading  of  tha  dally  observation 
tlaa. 

Lapse  rats  alteration  of  tha  BaxtBi»  dally  air  tnqpaxaturs. 

Lapse  rata  altarstloo  of  the  Binisn  dally  air  taaparatuxa. 

-2  -I 

Total  daily  ineoaing  solar  radiation,  MJn  day 
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VnHSH 

VQLOSH 

VSOIL 

HEl 

m(J,8) 

XXX(J,1) 

YEAR 

m(j.2) 

m<j,3) 

m<j,6) 

m{j,7) 

m<j,9) 

(2)  TSTM 

AB 

ACL(8} 

ALFB(IX) 

AO 

AFSH 

B 

BBB(J,I) 

BC 

BT  (8) 

BEP 


Victor  hoUlnc  ntm  anon  valuas  of  thaxaal.  dlffusivlty,  haat  conductivity, 
aadaiivity,  albodo  and  danaity. 

Vaetor  holdint  old  ano«  valuaa  of  thaiaal  diffuaivity,  haat  conductivity, 
aadaaivity,  alhado  and  danaity. 

Vaetor  holding  aoil  valuaa  of  thantal  diffuaivity,  haat  conductivity,  aoiiaaivity, 
alhado  and  danaity  (redundant). 

Uaad  in  tha  calculation  of  dapth  idian  converting  ’now  anow'  to  'old  anoo' ,  cm. 

Julian  data. 

Obaarvation  tlma. 

Yaar  of  tha  aiaulation  data. 

Mlniauai  air  taa^aratura,  *C,  tabla  1. 

Ralativa  huBoiidity,  X,  tabla  2 
Cloud  eovar,  0-1,  tabla  3. 

Hind  apaad,  ma  tabla  6. 

Fracipltation ,  ■■  watar. 

Hazimum  air  ti^araturo,  *C. 


Hugga-MBllar  abaorptlon  function. 

Coafficiant,  a,  dapandant  on  cloud  typa,  uaad  in  tha  calculation  of  CTF. 

-2  -1 

Thamal  diffuaivity  of  layer  IX,  cm  min 
Ataoapharic  alhado  for  Raylaigh  acattaring. 

PACTE»IEMP^,  ealem'*«in'^*C. 

-2  -1 

Boat  conductivity  of  aurfaca,  calca  min  *C. 

Y  Intaroopt  of  linear  oquatiem,  uaad  for  tabla  interpolation. 

Conatant  uaad  in  tha  calculation  of  watar  vapour  praaaura. 

Coafficiant,  b,  dapandant  cn  cloud  typo,  uaad  in  tha  calculation  of  CTF. 
Bottoai  boundary  thamal  IR  aaUaalvity. 
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BK 

BFBM 

BTERM 

CC 

COEl 

COE2 

COEF 

CLRd) 

CLOUD 

CP 

CTQU 

CTF 

DAY 

DCOEF 

DODT 

DECL 

DELI 

DEPTH(IY) 

DEPTB(4S0) 

DONRIR 

DPRM 

OFRMO 

Dimi 

DTEBM 

EA 

ELF 


Bottoa  ■uxfaea  (aoMtaie  ahapa,  fraction  (0. 0-1.0) 

Boat  cooduetlvltr  of  hottcm  boundary  layar. 

Enarsy  eontrlbutad  by  Inaolatlon  aftar  adjuataant  ualns  aurfaca  abaorptlvlty, 

,  *2  .  *1 
ealca  aln 

CLOUD  *  CLOUD. 

Coafflelant  uaad  In  tha  calculation  of  BTE8M,  DTERM,  EF  and  DTESH,  aat  according  to 
Rlchardaon'a  Huabar. 

Coafficlant  uaad  in  tha  calculation  of  BTERH,  DTESH,  EF  and  DTESH,  aat  according  to 
Rlcbaraon'a  Hu^ar. 

Fart  of  calculation  of  Intarfaca  valuaa. 

Coafficlant  dapandant  on  cloud  typo,  uaad  in  tha  calculation  of  incoming  longwave 
radiation. 

Cloud  eovar,  fraction  (0.1- 1.0). 

Specif ie  haat  of  dry  air  at  eonatant  praaaura. 

Air  ta^Mratura  in  dagraoa  kalvln,  uaad  in  calculation  of  avaporativa  heat  flux. 
Cloud  adjuatawnt  factor. 

Julian  day  uaad  in  aolving  inaolatlon. 

Part  of  calculation  of  intarfaca  valuaa. 
dT/dTg,  part  of  upper  boundary  calculation. 

Solar  docllnatlon  angle. 

Tima  atap  in  houra. 

Hatrix  bolding  layar  doptha. 

Depth  of  anow  ox  aoll  profile. 

-2  -1 

Incoa^ng  longwave  radiation,  calca  min 
-2  -1 

Haat  flux,  calem  min  ,  at  bottom  boundary  or  tamparatura  in  ratdeinaa  at  bottom 
boundary. 

Tamparatura  of  bottom  material,  *C,  uaad  whan  LFLUX-0. 

-2  -1 

Beat  flux  of  banaatb  bottom  material,  calem  min  ,  uaad  tdian  LFLUEY  not  equal  0. 
Energy  flux  duo  to  evaporation. 

Hatar  vapour  praaaura  of  tha  air. 

Latitude,  radlana. 
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EPSM 

ES 

ESAT(T) 

E(T) 

F2 

FACTA 

FACTO 

FACTE 

FACTB 

FK(IX) 

Ft*i(J,I) 

6 

GTERM 

HEADER 

HTERH 

lEFSWT 

lEOF 

III 

IMATL 

IHTRdX) 

IFRMT 

ITER 

ITIME 

IX 

IT 


blsslvtty  of  aurfaca  ■ataxtal. 

Hatot  v^oux  proaauxa  of  tlia  auxfaca. 

Satuxatod  vapoui  pxaaauxa  at  taa^axatuxa  T. 

Vapoux  pxaaauxa  at  taapaxatuxa  T. 

Paxt  of  tha  calculation  of  bottoai  boundaxy  valuaa. 

SIGMA*EFSM 

FACTIXIGMA*BK*BEP*TR**4 ,  uaad  In  bottoa  boundaxy  clculatlon  Mhan  thaxa  la  alrapaca 
banaath  tha  bottoai  layax. 

FACTE-SIGMA*BK*BEP,  uaad  In  bottoai  boundaxy  calculation  Mhan  thaxa  la  alrapaca 
hanaath  tha  hottoa  layax. 

FACTB- (1000.0 /PRESS) **0.286,  uaad  In  aolvlnc  convactlon  tarm  (HTERM) . 

Haat  conductivity  of  layax  IX,  calailn  ^cai  ^*K 

Slopa  of  llnaax  aquation,  uaad  fox  tabla  Intaxpolatlon. 

Accalaxatlon  dua  to  gravity. 

-2  -1 

Ground  haat  flux,  calca  min 

72  chaxactax  Input  vaxlabla  uaad  to  print  coananta  on  output  (not  utlllaad  In 

smo). 

-2  -1 

Enargy  loaa  ox  gain  dua  to  convactlon,  calcm  min 

Switch,  whan  -0  will  print  output  only  at  apaclflad  tlma,  If  not  -0  will  print 
output  at  avary  Itaratlon. 

Sat  from  0  to  1  whan  an  EOF  la  ancountarad,  uaad  to  tarmlnata  program. 

Count  of  numbar  of  Intarval  Itoratlcoa  In  calculation  of  uppar  boundary  valuaa. 
Backward  oountax  of  layaxa,  ataxtlng  with  tha  nuabar  of  layara. 

Baginning  aub- layax  dapth  nuadsar  fox  layax  nuoAiar  IX. 

Backward  countar  SET*BPRIIT .  Whan  -1  output  la  pxlntad. 

Itaratlon  eountax  uaad  In  flnlta  dltfaxanca  calculation  of  haat  flow  aquation. 
Backward  countax  Inltlallsa  aa  total  tlma  atapa  In  hour. 

Layax  nmdiax  ataxtlng  with  top  layax. 

Sub-layax  dapth  nimdiax. 
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JMAX 

KIO 

KSQ 

XRMPA 

XTIHFG 

LAYERS 

LAY 

LFLUXY 

Ul 

M 

NAX(J) 

HCLOUO 

NETRAD 

HIT 

RKOHT 

RTABL 

HX(IX> 

OUTCD 

PI 

PRESS 

PTYME 

QSATCT) 


Th*  total  nuBfaor  of  aub-Iayar*. 

Count  of  tba  maabar  of  Intarnai  Itacationa  of  TSIN,  In  ordar  to  caleulata  tha  daily 
■aans  of  tha  TSIM  ealeuatad  anarsT^budcat  and  taoq^ratura  varlablaa. 
von  Kataan'a  eonatant  aquarad. 

Air  tifiratura,  datraas  kalvin. 

Surfaca  taoparatura,  datxaaa  kalvin. 

Count  of  tba  profila  layara. 

Latituda  usad  in  aolvins  insolation . 

Input  bottoa  boundary  data  control  awitcb.  If  •O,  thara  is  no  boat  flux  through 
bottoB  of  aatarial,  if  nagativa  thara  ia  no  air  apaea  banaath  bottoa  nataxial,  if 
poaitiva  thara  ia  air  apaea  banaath  botton  natarial. 

Dmmr  variabla  to  road  lina  numbar  froai  input  fila. 

Saeant  of  aolar  zanitb  angla,  radians. 

Huid>or  of  input  tabla  valuas  uaod  in  tabla  intarpolation  oiodula. 

Cloud  typo  indax  nuiRioz  (1-9)  usad  in  solving  insolation,  IR  aolsalon. 

Rot  radiation,  (Un  ^day 

Htabar  of  itarationa  usad  in  boat  flow  calculations. 

Usad  with  KOHIRL  to  taralnats  SUCK). 

Tabla  nvadiar. 

HX(IX>-TBK(IX)/SFRQ(IX),  nuobar  of  sublayer  of  saeb  layer. 

Output  Manipulation  variabla. 

3.1S1S93 

Ataospharic  prassura,  sA,  used  in  solving  Insolation. 

Baginning  tisM  of  output<^otal  m^or  of  hours  Minus  21,  usad  in  print  output 
anduls  2. 

Saturated  apoclfic  hiaidity  at  taqporaturs  T. 

Spoclfic  biaidity  at  twapiratura  T. 


OCT) 


UP 

saoA 

RBOCCIXl 

RI 

U 

U(IX) 

SAZ 

SCF 

SFKQdX) 

SICF 

SIGMA 

SLOPE 

SMALLA 

SPEED 

STASH 

STORd.IY) 

STOR(2,IY) 

STORO.IY) 

SIGR(*,IY) 

STORC.IY) 

STCR(6,IY) 

STORCY.IY) 

SOC 

SUK 

SUUAC 

T 

TA 

TAC 

TAX 


Balislvlty  bmaatb  •irtpac*. 

Air  danalty,  caleulatad  aa  part  of  aanaibla  baat  calculation. 

FK(IX)/ALPB<IX),  eaL<ai'**K'^. 

Klehardaon'a  Htadiar,  Indaz  uaad  In  aolvlnc  BTEHM,  DTEBH,  BTER,  EF. 

Surfaca  banaath  airapaca  gaoaatrlc  thapa,  traction  (0. 0-1.0). 

BR(IX)~DELT/SFIIQ**2,  part  of  boat  flow  aquation. 

Solar  aaiauth,  radiana. 

Uaad  in  calculation  of  aanaibla  baat  flux,  dataninad  by  tba  Richardaon  nunbar. 

o  -1 

Vartical  grid  apaelng  in  cai  in  aacb  layar  IX,  ca  ■ 

Inaolatlon  adjuafant  dua  to  xanith  angla,  aurfaca  alopa  and  aurfaca  aapact  angla. 
Stafan-Boltaan  conataat,  8.12a'll. 

Surfaca  alopa,  dagraaa,  with  boriaontal  >0  dagraaa,  uaad  in  aolving  inaolation. 
Abaorbtivity  of  aurfaca  aatarial. 

Wind  apaad,  caaac 

Huwarical  atability  paraawtar,  SO.S  for  nwarical  atability. 

Eatiaata  aub-layar  taaparatura,  dagraa  rankina. 

FK,  baat  conductivity  of  aub-layar  lY,  calain  ^ca  ^*K 
HBOC,  FK/ALPB.  calca***x'^. 

Conatant  diaanaionlaaa. 

Initial  anow  tanparatura,  dagraa  rankina,  of  initial  anow  profila. 

Saaa  aa  8T0R(2.IY). 

SaM  aa  STCR(3,IY}. 

Count  of  dally  anargy  totala. 

Caleulatad  inaolatlon  valua. 

Surfaca  aslautb,  dagraaa,  with  aoutb  >0,  waat  >+90  and  aaat  >-90  dagraaa,  uaad  in 
aolving  inaolation. 

Saaw  aa  IDC. 

Air  taaqiaratura.  dagraaa  rankina. 

Air  taaparatura,  *C. 

Air  taaparatura,  ’K. 


TAL 


Conatant  uaad  in  tha  ealoulatloo  of  HATER. 


TB 

ID 

TBC 

TEMi 

TFRQ 

THK(IX) 

TIME 

TIMER 

TITU 

lOITIM 

TFRNT 

IR 

TSK 

TIME 

HATER 

MET 

XCXIJ.l) 

XXX(J,2) 

XXX(J,3) 

XXX(J,4) 

XXX(J.S) 

TTYCJ,*) 

Z 

ZA 

ZASB 

ZZA 
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Thaiaal  eanductlvlt;  of  botton  Mtarlal,  calea  ’C  aln 
Daw  point  taapatatuia  usad  In  tha  calculation  of  HATER. 

Surfaca  taaparatura  of  aatacial.  dagraaa  ranklna. 

Bottoa  lapax  taaparatura  of  aatarlal  In  dapraoB  ranklna. 

Tlaa  atap  In  alnutaa  usad  In  aolvlng  hast  flow. 

Layar  thlcknaas.  In  ca  of  layax  IX. 

Tlaa  In  houxa  tihlch  aatarlal  taaparaturaa  ara  aatlawtad. 

Sun 'a  hour  an(la.  cadlana. 

Surfaca  taaparatura ,  *C . 

Total  nvafaar  of  24  hour  rapatltlona  uaad  In  aolvlng  hast  flow. 

Output  tlaa  print  fraquancy,  ainutaa. 

Taaparatura  of  alrapaca  banaathbottoa  aatarlal. 

Matarlal  aub- layar  taaparatura,  *K. 

Tlaa  In  boura,  uaad  In  Inaolatlon  calculation. 

Tha  aaount  of  praelpltabla  watar,  aa,  uaad  in  aolvlng  Inaolatlon. 

Nolatura  contant  of  aurtaca  aatarlal. 

Tlaa  in  hours  for  tabla  1  (air  taagparstura). 

Tlaa  In  hours  for  tabla  2  (ralatlva  bualdity). 

Tlaa  In  hours  foe  tabla  3  (aaount  of  cloud  covar) . 

Tlaa  In  hours  for  tabla  4  (solar  Insolation). 

Daptb,  ca,  for  table  S  (taaparatura  profile)  Initial  taaparatura  profile, 

-2  -1 

Insolation,  calea  ain  ,  If  0.0  at  1200  boura,  insolation  values  will  be 
calculated,  table  4. 

Solar  sonlth  angle. 

Shelter  height,  ca. 

Shaltar  height,  ca. 

Surfaca  taaparatura  of  aatarlal,  degree  ranklna. 


ZZB 


Bottoa  layar  taaparatura  of  aatarlal  in  dagraa  ranklna. 


ATERH 


-2  *1 

Ataaaphaxie  IR  aalaalon,  caica  min 

ATF(l)  Fart  of  tha  root-flndlnt  alsorlthB. 

ATF(2)  Part  of  tha  root'tlndlnt  alRorlttaa. 

BINT  Part  of  tha  root-tindinc  alcotitha. 

CBO  Dlaanaloolaaa  haat  or  aolatura  tranafar  coafflclant  appllcabla  to  tha  top  of  a 

danaa  canopy. 

CBH  Dlaanalonlaaa  haat  or  anlatura  tranafar  coafflclant  applicable  to  tha  top  of  a 

danaa  canopy. 

CBS  Baat  tranafar  coafflclant  aquation  Intarpolatad  batwaan  ground  with  no  covac,  CBO 

and  coaplata  cover,  CBB. 

CP  Specific  haat  of  dry  air  at  conatant  praaaura. 

DELTMF  Incraaantal  tnapiratura  for  tha  root-finding  algorltha. 

DTBETA  DTHETA>(TA-TF)*FACTB/ZA.  uaed  In  the  calculation  of  tha  Rlchardaon  Mianbar. 

DU  Uaad  In  tha  calculation  of  tha  Rlchardaon  Ruabar. 


E(T) 

EF 

EFF 

ELG 

EPl 

EPF 

EPSB 

ESAT(T) 

EX 


Vapour  praaaura  at  tanparature  T. 

Latent  haat  loaa  to  the  atooaphara  for  tha  foliage. 

Effective  tanparatura. 

Latent  heat  loaa  to  tha  ataoaphare  at  the  ground  surface. 

Fart  of  tha  calculation  of  tha  energy  budget  for  foliage . 

Balsslvlty  of  foliage. 

Egdsslvlty  of  ground  (snowpack)  surface. 

Saturated  vapour  praaaura  at  taovaratura  T. 

Coefficient  uaad  In  tha  calculation  of  BTERM,  DTERH,  BTER  and  EF,  sat  according  to 
Richardson's  Hiabar. 


FACTS  Potential  taaparatura. 

FEB(l)  Part  of  tha  root-finding  algorlthB. 

Part  of  tha  root-finding  algorltba. 


FEB(2> 
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FEHB  En*xgy-bud(«t  of  tho  ground  and  folios*  surfocos  rospoctivaly. 

FIC(IX)  Hoot  conductivity  of  (round  (snowpack). 

FOLA  Abaocptlvlty  of  th*  follas*  layer  (l-all>*do>. 

-2  -1 

FOLQB  Follat*  thamal  IS  aalttane*.  calca  nln 

G  Aeealaratlon  dua  to  (ravtty. 

-2  -1 

GBMD6B  Ground  thamal  IB  anlttanea,  calca  Bin 

HFOL  Folios*  halsht,  ca. 

BSF  Convactlv*  anorsy  flux  acroa*  th*  follas*  aurfac*. 

BSG  Convactlv*  anarsy  flux  acroas  th*  (tound  (aaowpack)  aurfac*. 

HTER  Part  of  th*  calculation  of  BSF. 

OAF  Spaclfle  huaddity  of  air  within  th*  follas*  layer. 

OF  Specific  hunldlty  of  th*  air  at  follas*  halgbt. 

06  Spaclfle  hualdlty  of  th*  (anowpack)  aurfac*. 

QSAT(T)  Saturated  apaclfle  hmldlty  at  taoswratur*  I. 

NDEX  Count  of  th*  nu^ar  of  tinea  th*  *D*rsy~buds*t  for  the  ground  or  tb*  foliage  was 

calculatad. 


SA 

SC 

BDF 

RZFRAD 

BBOAF 

BBQAG 

KI 

BID 

BLU 

BS 


Atnoapharlc  raalatanc*. 

Canopy  raalatanc*  to  water  vapour  dlffualon 

Fraction  of  potential  evaporation  rat*  freo  follag*. 
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Surface  related  radiation,  calcm  nln 
Danalty  of  air  near  ground  aurfac*. 

Danalty  of  air  near  th*  ground  aurfac*. 

Rlehardaon  Rtabar. 

Downward  longwav*  flux  ovar  ground  (snowpack)  surface. 

Upward  dlractad  longwave  flux  radlatad  fren  ground  (snowpack)  aurfac*. 
Stoawtal  raslatane*. 


SPRQdX)  Dlatane*  of  first  grid  point  below  th*  aurfac*. 


S6 


Ineonlns  solar  radiation  raachlns  the  ground  (snowpack). 


SHRH 

8I(V 

SIGMA 

SLOFBl 

SMALLA 

SOL 

STATE 

SURFOB 

T1 

TA 

TEFF 

TEFFR 

TEML 

TF 

TFO 

TF* 

TG* 

THETAV 

OA 

OAF 

NET 

ELI 

xunw 

EO 

ZA 

ZD8F 
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ABooBt  of  ahortHav*  radiation  abaorbad  by  toltaaa. 

Follata  eoaar  traetlan. 

Stafai-Boltanann  eonatant. 

Fart  of  tha  root-fladlac  aigorithB. 

Ground  (anoupaok)  aurfaoa  abaorptivlty. 

Caleulatad  aolar  inaolation  raaehlnt  tba  top  of  tha  follaga  layar. 

Arbitrary  nultlpUar  (STATE><))  of  RS  uaad  to  account  for  sancaaanca,  atraaa  ate. 

-2  -1 

Surfaca  (foUaga  *  grotmd  aurfaca)  tbaraal  IR  anittanca,  calca  aln 

Tanparatura  at  flrat  grid  point  balow  aurfaca,  dlatanca  SSFQ(IX),  in  TSIH  boat 
tranafar  algorltbai. 

Air  tanparatura. 

Moan  affactlva  tanparatura. 

TEFF  with  raflactloo  Incorporatad. 

Graiad  (anowpack)  aurfaca  taaqparatura. 

Follaga  tanparatura. 

Fart  of  tha  root-finding  algorltiB. 

Fart  of  tha  caloulatloo  of  FERB. 

Fart  of  tha  oaloulatlon  of  FERB. 

Uaad  in  tha  calculation  of  Riebardaon'a  Hindiar. 

Wind  apoad. 

Wind  apaad  of  air  In  tha  follaga  layar. 

Nolatura  contant  of  groimd  (anowpack)  aurfaca. 

L~f(T^),  uaad  In  tha  oalculatloo  of  tha  latent  baat  of  evaporation,  a  function  of 
air  taa^raturo. 

Fart  of  tha  foliage  anargy-budgat  calculation. 

Roughnaaa  langth. 

Inatnaant  balgbt  abova  ground  (anowpack). 

Zero  dlaplaeaaMBt  balgbt. 
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